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An  exact  theory  regarding  solitary  internal  gravity  waves  in  stratified 
fluids  is  presented.  Two-dimensional,  inviscid,  incompressible  flows  confined 
between  plane  horizontal  rigid  boundaries  are  considered.  Variational 
techniques  are  used  to  demonstrate  that  the  Euler  equations  possess  solutions 
that  represent  progressing  waves  of  permanent  form.  These  are  analogous  to 
the  surface,  solitary  waves  so  easily  generated  in  a  flume.  Periodic  wave 
trains  of  permanent  form,  the  analogue  of  the  classical  cnoidal  waves,  are 
also  found.  Moreover,  internal  solitary-wave  solutions  are  shown  to  arise  as 
the  limit  of  cnoidal  wave  trains  as  the  period  length  grows  unboundedly. 
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SIGNIFICANCE  AND  EXPLANATION 


The  notion  of  a  solitary  wave  arose  in  the  last  century  and  provoked  some 
controversy  as  to  the  existence  of  such  a  phenomenon.  After  more  than  half  a 
century  of  only  sporadic  interest  in  these  waves,  the  last  twenty  years  have 
seen  an  upsurgence  of  scientific  work  on  the  subject.  The  interest  stems  from 
the  realization  that  these  special  wave  forms  play  a  significant  role  in  the 
evolution  of  general  classes  of  disturbances.  This  property  of  solitary  waves 
is  observed  on  a  particularly  grand  scale  in  the  earth's  oceans.  There, 
disturbances  that  appear  to  be  internal  solitary  waves  with  crests  spanning 
hundreds  of  kilometers  have  been  recorded.  Various  approximate  models  have 
been  used  to  analyze  these  motions.  Nevertheless,  we  are  far  from  a  complete 
mathematical  treatment  of  the  phenomena  just  described. 

In  this  report  a  model  physical  problem  is  studied  in  a  mathematically 
exact  formulation.  We  consider  two-dimensional  flows  confined  between  rigid 
horizontal  boundaries  and  show  that  there  are  solutions  of  the  Euler  equations 
representing  internal  solitary  gravity  waves.  Our  theory,  which  is  not 
restricted  to  small  amplitudes  predicts  both  waves  of  elevation  and 
depression,  depending  on  the  ambient  density  distribution  and  the  velocity 
distribution  at  infinity.  Just  as  for  the  classical  surface  solitary  waves, 
these  waves  are  single-crested,  symmetric,  and  decay  exponentially  away  from 
the  crest.  They  thus  represent  disturbances  of  essentially  finite  extent. 
These  qualitative  features  are  established  using  variational  techniques 
combined  with  the  theory  of  rearrangements.  Moreover,  it  is  shown  that  the 
single-crested  wave  arises  as  a  limit  of  periodic  wave  trains  of  increasing 
wave  lengths. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


FINITE-AMPLITUDE  STEADY  WAVES 
IN  STRATIFIED  FLUIDS 

J.  L.  Bona*'1'2,  D.  K.  Boset+'3  and  R.  E.  L.  Turner*'1'4 

1 .  INTRODUCTION 

The  mathematical  problem  to  be  analysed  in  this  paper  arises  in  the  study 
of  steady,  two-dimensional  wave  motion  in  heterogeneous,  inviscid, 
incompressible  fluids  confined  between  two  rigid  horizontal  planes.  The 
existence  and  various  properties  of  f inite-amplitude  periodic  wave  trains,  and 
single-crested  waves,  of  permanent  form  will  be  established.  These  internal 
waves  are  analogous,  respectively,  to  cnoidal  and  solitary  surface  waves,  and 
will  be  given  these  names  in  what  follows. 

The  first  rigorous  existence  theory  for  internal  waves  appears  to  have 
been  presented  by  Kotschin  (1928).  He  considered  a  system  comprised  of  two 
homogeneous  layers  of  fluid  of  differing  density,  the  lighter  resting  on  the 
heavier,  and  the  whole  contained  between  two  fixed  horizontal  plane 
boundaries.  Kotschin  used  a  major ant  method  to  obtain  small -amplitude 
cnoidal-wave  solutions  of  the  Euler  equations  posed  by  this  two-fluid 
system.  The  first  rigorous  theory  to  include  both  cnoidal  and  solitary 


t 

Department  of  Mathematics,  University  of  Chicago  and  the  Mathematics 
Research  Center,  University  of  Wisconsin-Madison,  U.S.A. 
ft 

Department  of  Mathematics,  Brighton  Polytechnic,  U.K. 

* 

Department  of  Mathematics  and  the  Mathematics  Research  Center,  University  of 
Wisconsin-Madison,  U.S.A. 

- 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 ; 

2 

the  National  Science  Foundation  (U.S.A.)  under  Grant  No.  MCS-8002327; 

3 

the  Science  Research  Council,  U.K.; 

4 

the  National  Science  Foundation  (U.S.A.)  under  Grant  No.  MCS-7904426. 


internal  waves  was  given  by  Ter-Krikorov  (1963).  He  considered  continuously 
stratified  fluids  with  a  rigid  bottom  boundary  and  either  a  fixed  or  a  free 
top  boundary.  Ter-Krikorov  relied  on  Yih's  (1960)  version  of  long's  (1953) 
equation  and  on  the  approximate  description  of  the  waves  in  question  that  are 
valid  for  small  amplitudes  and  long  wavelengths.  Casting  the  equation 
satisfied  by  the  difference  between  the  exact  and  approximate  solutions  into 
operator  form,  he  applied  the  contraction-mapping  theorem  to  this  operator, 
considered  as  a  mapping  of  a  small  ball  centered  at  the  zero  function  in  an 
appropriate  function  space.  In  this  manner,  he  was  able  to  deal  with  both 
cnoidal  and  solitary  internal  waves,  but  only  of  very  small  amplitude.  His 
work  thus  compares  with  that  of  Lavrentief  (1943,  1947),  Friedrichs  and  Hyers 
(1954)  and  Littman  (1957)  concerning  surface  waves  of  permanent  form  in  liquid 
of  finite  depth.  Zeidler  (1971),  using  complex  function  theory,  reduced  the 
steady  internal-wave  problem  in  the  above-described  two-fluid  system  to  a 
situation  in  which  he  could  demonstrate  that  a  bifurcation  at  a  simple 
eigenvalue  occurs.  The  bifurcating  branch  thus  gave  small-amplitude  cnoidal 
waves,  the  general  result  being  closely  allied  in  spirit  to  that  of 
Kotschin.  Going  beyond  Kotschin's  theory,  Zeidler  was  able  to  take  account  of 
capillarity  in  his  theory.  Benjamin  (1973)  gave  an  exact  treatment  of  cnoidal 
internal  waves  whose  amplitudes  need  not  be  small.  He  made  use  of  positive- 
operator  methods,  as  pioneered  by  Krasnosel 'skii,  together  with  the  Leray- 
Schauder  degree  theory.  His  method  did  not  yield  a  theory  for  solitary  waves, 
though  some  mathematical  evidence  for  their  existence  was  presented. 

The  present  treatment  of  the  problem  is  very  much  influenced  by 
Benjamin's  (1973)  wide-ranging  article,  though  our  methods  are  different. 
Benjamin's  analysis  has  the  speed  of  propagation  as  a  fixed  parameter  whilst 
the  energy  possessed  by  the  wave  is  left  free.  This  problem  is  covered  here. 


as  well  as  the  complimentary  specification  in  which  the  energy  is  fixed  and 
the  speed  of  propagation  is  left  free.  As  explained  in  section  7,  this  latter 
case  is  probably  more  relevant  to  typical  situations  in  which  such  waves  are 
generated.  Our  variational  methods,  when  combined  with  the  use  of 
rearrangement  inequalities,  are  effective  for  a  broad  range  of  non- 
linearities.  This  combination  of  techniques  has  been  fruitfully  used  in  other 
hydrodynamic  problems  (see  e.g.  Friedrichs  1934;  Garabedian  1965;  Fraenkel  and 
Berger  1973)  and  our  treatment  of  constrained  variational  problems  has  several 
points  in  common  with  the  last  article  cited.  The  use  of  symmetrization  in 
the  problem  of  fixed  speed  appears  to  be  new. 

In  an  outgrowth  of  the  present  paper,  Turner  (1982)  has  used  an  alternate 
variational  principle  to  establish  the  existence  of  cnoidal  and  solitary  waves 
in  a  fluid  with  a  rapidly  varying  heterogeneity,  the  two-fluid  system 
mentioned  earlier  being  a  special  case  of  his  theory. 

The  plan  of  the  paper  is  as  follows.  Section  2  is  devoted  to  a  general 
description  of  the  mathematical  problem  that  presents  itself  in  considering 
internal  waves  of  permanent  form.  In  section  3,  existence  of  cnoidal  internal 
waves  is  established,  both  for  the  case  where  the  speed  of  propagation  is 
given  and  for  the  situation  in  which  the  total  energy  is  fixed,  h  priori 
bounds  satisfied  by  these  periodic  solutions  are  deduced  in  section  4.  These 
bounds  are  independent  of  the  period  length,  and  include  a  result  of 
exponential  decay  from  crest  to  trough.  In  section  5  internal  solitary  waves 
are  shown  to  arise  as  the  limiting  forms  of  internal  cnoidal  waves  as  the 
period  length  becomes  indefinitely  large.  Also  discussed  is  a  proposition 
implying  the  absence  of  closed  streamlines.  This  technical  point  is  important 
for  interpreting  the  solutions  obtained  in  our  analysis  as  realizable  wave 
motion.  A  broad  class  of  concentration  profiles  is  examined  in  section  6,  and 


shown  to  yield  equations  that  fall  within  the  confines  of  our  theory.  The 


relationship  of  our  theory  to  other  studies  is  considered  in  section  7.  Field 
and  laboratory  evidence  concerning  internal  waves  of  permanent  form  is  briefly 
reviewed.  Various  implications  and  drawbacks  of  the  theory  developed  herein 
are  presented,  along  with  suggestions  of  avenues  for  further  investigation. 


2.  THE  GOVERNING  EQUATIONS 

To  understand  certain  aspects  of  the  theory  to  be  developed  presently, 
the  derivation  and  import  of  the  equations  that  will  hold  our  attention  must 
be  kept  in  mind.  It  is  therefore  useful  to  review  some  points  arising  in  the 
formulation  of  the  mathematical  problem  to  be  analysed.  The  description  given 
here  follows  closely  Benjamin's  (1973)  account,  but  the  gist  of  it  goes  back 
much  further  (cf.  Benjamin  1966,  Dubreil-Jacotin  1935,  U>ng  1963  and 
especially  Yih  1958,  1960). 

The  idealized  physical  situation  to  be  modeled  is  a  heterogeneous 
incompressible  fluid  moving  between  horizontal  plane  boundaries.  Interest 
will  be  focused  on  two-dimensional  flows,  which  will  be  assumed  to  be  inviscid 
and  non-diffusive.  These  flows  will  be  described  in  a  two-dimensional 
Cartesian  coordinate  frame,  in  which  the  x-axis  lies  along  the  bottom  flow 
boundary  and  height  y  above  this  boundary  comprises  the  second  coordinate. 
The  distance  between  the  boundaries  is  taken  as  the  unit  of  length,  so  the 
lower  and  upper  boundaries  are  represented  by  {(x,0)  s  x  e  r}  and 
{ (x, 1 )  :  x  e  R},  respectively. 

There  is  postulated  a  primary  or  base  flow  in  which  the  fluid  motion  is 
everywhere  horizontal  with  velocity  U(y)  at  height  y  above  the  lower 
boundary,  0  <  y  <  1.  It  is  also  assumed  that,  in  the  primary  state,  the 
fluid  density  p  is  a  decreasing  function  of  y,  for  0  <  y  <  1.  In  fact, 
it  will  be  assumed  that  P  is  a  continuously  differentiable  function  of  y, 
and  that  P'(y)  <  0,  for  y  in  [0,1]. 

We  search  for  waves  of  permanent  form  whose  velocity  of  propagation 
downstream  (in  the  direction  of  increasing  x)  is  c.  It  is  convenient  then 
to  presume  that  the  coordinate  system  is  also  moving  downstream  at  speed  c, 
thus  rendering  stationary  the  waveform  in  question.  In  the  moving  frame  of 
reference,  the  primary  fluid  velocity  is  W(y)  =  U(y)  -  c. 
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Let  g  =  (u,v)  denote  the  velocity  field  of  a  steady  flow  of  interest. 

The  flow  is  assumed  to  be  incompressible  and  non-dif fusive.  The  former 

assumption  means  g  is  divergence  free,  whilst  the  latter  amounts  to  the 

assertion  that  the  derivative  of  the  density  following  the  flow  is  zero.  Thus 

V»g  =  0  and,  since  the  flow  is  steady,  g*Vp  =  0.  It  follows  that 
1  /2 

V»(p  g)  =  0,  so  there  is  a  pseudo-stream-function  ifi  such  that 

1/2 

(♦„,  -  ♦„>  =  P  3  ,  (2.1) 

y  x 

and  it  may  be  supposed,  by  way  of  normalization,  that  ip(x,0)  =  0.  It  follows 
readily  that  P  =  P(^)  is  a  function  of  alone.  Note  that  in  the  special 
case  of  a  parallel  base  flow  g  =  (W(y),0),  the  associated  normalized  pseudo¬ 
stream-function  is 

¥{y)  -  !yQ  PV2(z)W(z)dz  ,  (2.2) 

with  P  being  the  specified  function  of  height  above  the  bottom  boundary. 

The  Euler  equations  for  steady  flow  are 

P(g*7)g  =  -?p  -  gpk  ,  (2.3) 

where  p  is  the  pressure,  g  the  gravity  constant,  and  k  the  unit  vector 
(0,1).  Keeping  in  mind  (2.1)  and  the  dependence  of  P  solely  on  «|»,  and 
introducing  n  =  Ai|>#  we  easily  derive  from  (2.3)  the  relations, 

^♦x  =  3X<P  +  \  Pl^l2)  > 

and  (2.4) 

=  9w(p  +  4  Ptsil2)  +  9P  • 

y  y  ^ 

1  .  .2 

Let  H  =  p  +  —  P|g|  +  gPy  be  the  total  head  associated  to  our  steady  flow. 
Now  H  is  also  a  function  of  ip  alone,  as  follows  immediately  from  (2.4)  and 
the  fact  that  P  =  P(^).  If  H  is  differentiated  with  respect  to  <|»  using 
the  chain  rule,  and  the  result  simplified  by  employing  the  relations  (2.4), 
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there  appears 


+  H'(^)  =  gyP '  ( ty)  ,  (2.5) 

where  the  '  denotes  differentiation  with  respect  to  ||>.  This  is  Yih's 
version  of  Long's  equation  (cf.  Dubreil-Jacotin  1935,  Long  1953  and  Yih  1958). 
The  auxiliary  stipulations  are  the  kinematical  conditions, 

>|>(x,0)  =  *(0)  =  0  and  4»(x,1)  =  'i'(l)  ,  (2.6) 

for  x  e  R,  and  the  asymptotic  conditions, 

♦  (x,y)  ♦  'P(y),  as  |x|  •  ,  (2.7) 

for  0  <  y  <  1 .  Thus  ^  represents  a  flow  that  is  connected  to  the  primary 
state  at  infinity,  and  for  which  p  and  H  are  constant  along  its  pseudo¬ 
streamlines.  The  situation  envisaged  is  illustrated  in  figure  1. 

In  (2.5)  neither  P  nor  H  is  immediately  known  as  a  function  of  i|». 
Because  of  (2.7),  they  may  be  determined  in  principle  from  the  primary  flow. 
For  example,  if  U(y)  =  d,  a  constant,  and  c  =  d-c,  then  (2.2)  yields 

¥(y)  m  c  /q  p,/,2(z)dz  . 

Since  P  >  0,  f  is  strictly  increasing  and  so  may  be  inverted.  Let 
-1 

Y  =  "  .  Then  the  density  p  is  expressed  as  a  function  of  the  pseudo¬ 

stream-function  value  Z  by 

P(Z)  =  P ( Y ( Z )  )  . 

Since  P  is  constant  along  pseudo-stream-lines,  its  value  at  a  point 
(x0,y0),  in  the  flow  corresponding  to  ij>,  is  determined  by  tracing  the 
.'Udo-steamline  with  value  Wxg'yg)  to  infinity,  where,  because  of  (2.7), 
p  is  determined  by  the  last  displayed  relation.  Note  that  once  P(40  is 
known,  H  is  determined  from  (2.1)  and  the  hydrostatic  pressure  law  dp/dy  = 
-gP  for  the  primary  flow  . 

Define  the  perturbed  pseudo-stream-function  <j>  =  <t>(x,y)  by  the  formula, 

iMx,y)  =  'f(y)  +  j$(x,y)  ,  (2.8) 
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where  j  is  a  normalizing  constant.  The  boundary-value  problem  expressed  in 
(2.5,  2.6,  2.7)  reduces  to  the  nonlinear  eigenvalue  problem, 

+  h(y,4>)  =  Xf  ( y ,  4> )  (2.9a) 

with  the  supplementary  conditions 

${x,0)  =  $(x,1)  =  0,  for  x  6  R 

and  \  (2.9b) 

4>(x,y)  +  0,  as  |x|  “  ,  J 

for  all  y  in  [0,1].  The  eigenvalue  X  is  proportional  to  gc  ,  where 
c  is  a  velocity  scale  for  the  primary  flow  in  the  travelling  coordinates 

introduced  earlier  (e.g.  c  •  sup  W(y)).  Just  as  for  p(i]0  and  H(<|>), 

0<y<1 

the  functions  h(y,$)  and  f(y,$)  are  in  general  a  little  complicated  to 
determine  from  a  given  base  state,  and  this  aspect  is  ignored  here.  It  is 
dealt  with  in  detail  in  section  G,  where  a  set  of  examples  will  be  exhibited 
and  shown  to  fall  within  the  class  of  equations  to  which  our  theory  is 
applicable. 

The  remainder  of  this  section  is  devoted  to  various  aspects  of  the 
boundary-value  problem  (2.9)  and  its  set  of  solutions.  First,  by  the  term 
solitary  wave  we  will  understand  a  solution  $  of  (2.9)  which  is  even  in 
x,  monotone  for  x  >  0,  and  rapidly  convergent  to  zero  as  |x|  ♦  00 .  The 
associated  pseudo-stream-function  =  ’i'  +  j<J>  represents  a  steady  flow,  as 
described  previously,  which  is  a  single-crested  wave  of  translation,  symmetric 
about  its  crest,  and  decaying  rapidly  to  the  base  flow  away  from  its  crest, 
thus  being  in  essence  a  disturbance  of  finite  extent.  The  description  of  such 
flows  as  solitary  waves  is  provisional,  however.  In  order  that  such  waves 
fall  within  the  presently  accepted  conception  of  solitary  waves  (cf.  Scott 
Russell  1844  and  Miura  1976)  their  stability  must  be  ascertained. 
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Similarly,  solutions  $  of  (2.9a)  which  vanish  at  y  =  0  and  y  =  1, 
and  which  are  even  periodic  functions  of  x,  say  of  period  2k,  and  which 
are  monotone  on  tO,k]  will  be  referred  to  as  internal  cnoidal  waves.  The 
name  is  chosen  in  analogy  with  the  classical  two-dimensional  waves  on  the 
surface  of  a  liquid  (cf.  Kortweg  and  de  Vries  1895).  For  it  is  imagined  that 
the  flow  corresponding  to  such  a  $  is  a  periodic  wavetrain,  propagating 
steadily  and  without  change  of  form,  and  single-crested  within  its  fundamental 
period.  The  reader  is  cautioned  that  this  interpretation  cannot  be  strictly 
justified.  For,  according  to  the  derivation  of  Yih's  equation  (2.5),  unless 
the  pseudo-streamlines  connect  asymptotically  to  those  of  the  base  flow  (i.e. 

♦  0  as  | x |  ®),  the  functions  p  and  H  cannot  be  evaluated  by  passing 

far  upstream  as  suggested  earlier  (though  we  will  see  that  for  large  periods, 
the  flow  in  the  trough  agrees  essentially  with  the  primary  flow). 

Consequently,  whilst  a  periodic  solution  $  of  (2.9a)  does  indeed  define  a 
solution  of  the  Euler  equations  (2.3),  it  may  not  correspond  to  a  perturbation 
of  the  given  base  flow.  Its  direct  hydrodynamical  interpretation  thus  becomes 
clouded.  A  more  subtle  interpretation  of  the  generation  and  role  of  such 
periodic  internal  wavetrains  has  been  proposed  by  Benjamin  (1966),  and  this 
possibility  receives  some  attention  in  section  7. 

The  issue  of  whether  or  not  the  pseudo-streamlines  each  connect  to  a 
definite  value  at  infinity,  as  in  figure  1,  arises  even  for  the  solitary-wave 
solutions  of  (2.9).  For  solutions  of  (2.9a)  which  satisfy  the  boundary 
condition  (2.9b)  might  conceivably  take  values  strictly  outside  the  range 
[0,T(1)].  If  this  occurs,  there  would  necessarily  be  closed  streamlines  in 
the  finite  region  of  the  flow  domain,  representing  a  motion  with  an  eddy  of 
fluid  not  present  in  the  base  flow.  Indeed,  closed  eddies,  often  called 
"rotors'  do  appear  in  practice,  associated  with  larger  amplitude  waves,  as  in 


9 


9 
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the  experimental  observations  of  Davis  and  Acrivos  (1967).  It  is  plausible 
that  the  fluid  encompassed  within  a  rotor  is  nearly  homogeneous,  reflecting 
the  mixing  that  occurs  in  such  a  structure.  One  might  therefore  be  tempted  to 
model  this  situation  by  extending  the  functions  P  and  H  outside  the  domain 
[0,f(1)]  with  the  specification,  P(x)  >»  P(0),  for  <  0,  and 
P(x)  =  p(Y(1)),  for  >  I'd),  and  similarly  for  H  (cf.  Benjamin  1971). 

In  view  of  the  serious  limitations  of  perfect-fluid  theory  in  describing  this 
sort  of  flow  feature,  some  circumspection  regarding  such  eddies  3eems 
warranted.  Accordingly,  in  section  5  a  theorem  is  proved  which  may  be 
interpreted  as  stating  that  closed  eddies  do  not  occur  in  a  range  of  the  flows 
whose  existence  is  established  by  our  theory. 

The  final  point  respective  of  (2.9)  concerns  the  possibility  of  trivial 
solutions  of  (2.9a).  The  function  $  =  0  is  a  solution  of  (2.9)  that 
corresponds  to  the  undisturbed  base  flow.  As  it  happens,  there  is  generally 
another  trivial  solution  called  the  conjugate  flow  (cf.  Benjamin  1971).  This 
is  the  unique,  non-zero,  non-negative  x-independent  solution  of  (2.9a)  that 
respects  the  first  boundary  condition  in  (2.9b).  The  conjugate  flow  plays  an 
important  role  in  the  interpretation  of  various  wave  phenomena.  However,  as 
the  reader  will  readily  appreciate,  these  two  trivial  solutions,  the  base  flow 
and  its  conjugate  flow,  potentially  complicate  a  search  for  non-trivially  x- 
dependent  solutions  of  (2.9)  that  are  of  especial  interest  here. 


3.  EXISTENCE  OF  PERIODIC  WAVES 


In  this  section  we  prove  the  existence  of  solutions  (A,$)  of  the 
boundary  value  problem, 


-A«ji(x,y)  +  h(y,<j>(x,y) )  =  Af  ( y ,  $(x,y ) ) 
$(x,0)  =  <Mx,1)  =  0  , 


(P) (3.1) 


which  are  periodic  in  x,  where  ft  =  {(x,y)|-°°  <x<®,  0<y<  l}.  In  the 
problem  (3.1)  we  assume  that  f  and  h  are  defined  for  all  values  of  As 

we  ultimately  show  that  solutions  4>  of  (3.1)  determine  pseudo-stream- 
functions  ♦  whose  values  lie  in  the  physical  range  [0,  ¥(1)]  (cf. 
corollary  4.5)  this  embodies  no  essential  restriction.  The  basic  assumptions 
on  f  and  h  are  contained  in 


Condition  A. 


AI 


The  function  f  has  the  form. 


f(y,t) 


tfQ(y)  +  f1(y/t),  for  0<y<1,  t>0 


-f(y,  -t),  for  0<y<1,t<0, 


where  fQ,  f^  are  HOlder  continuous  on  bounded  sets;  fp  >  0  on 
0  <  y  <  1;  and  fj(y,t)  *  o(t),  uniformly  for  0  <  y  *  1,  as  t  +  0. 

Further,  there  exist  constants  m>1,n>1,  a>0  and  d  >  0  such  that  for 
t  >  0 

«tm  <  f  (y,t)  <  d( 1  +  tn)  . 

The  function  h  has  the  form. 


All 


h(y,t) 


C  thQ(y)  +  h1<y»t>'  for  o  <  y  <  i,  t>o, 

1  -h (y,  -  t),  for  0  <  y  <  1,  t<0. 
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where  hg  and  are  Hfllder  continuous  on  bounded  sets,  ^  is  continuous 

3h  2 

and  >  e  for  some  constant  e  >  -»  •  Further  there  aro  constants 
a  >  0  and  o*  such  that  for  t  >  0, 

o't”  <  h^y.t)  <  at”, 

and 

|h1 (y,t) |  <  d{ 1  +  tn), 

where  m,  n,  and  d  are  the  constants  appearing  in  AI.  Note  that  the  first 
inequality  above  implies  h.j(y,t)  *  o(t),  uniformly  for  0  <  y  <  1,  as 
t  ♦  0. 


The  final  assumption  and  the  subsequent  variational  principles  involve 
the  functions 


F(y,t)  -2  /*  f(y,s)ds,  F^y^t)  “2  /*  f^y^ds 
H(y,t)  «  2  /*  h(y,s)ds,  H^y^)  -  2  /Jj  h^y^sjds  . 


(3.2) 


The  last  hypothesis  is, 

AIII  For  each  X  >  a  (a  and  o  the  constants  in  AI  and  All) 
there  is  a  0  <  1  such  that 

XFt  (y,t)  -  H^yft)  <  ©(Xf^y^Jt  -  h^yjtjt)  . 


While  the  assumptions  are  rather  technical,  they  are  tailored  with  two 
ends  in  view.  First  they  are  weak  enough  to  be  readily  verified  in  a  range  of 
applications.  Further,  they  are  strong  enough  to  mitigate  certain  technical 
difficulties  associated  with  our  proofs. 

The  problem  before  us  is  to  find  solutions  of  problem  (P)  which  are 
periodic  in  x  with  period  2k.  We  consider  two  formulations  of  the 
problem.  The  first  is  a  constrained  problem  (PC): 


3- 


solve  (P) 


subject  to  /J  /*  [|V$|2  +  H  ( y ,  ♦ )  j dxdy  =  R2,  \  (PC)(3.3> 

where  R  is  a  given  positive  constant. 


The  second  is  a  free  problem  (PF): 


where 


solve  (P) 

X  is  a  given  constant. 


} 


(PF) (3.4) 


The  problem  (PC)  corresponds  to  specifying  the  "energy"  of  a  wave,  while 
(PF)  corresponds  to  specifying  its  velocity. 

The  analysis  of  both  problems  (PC)  and  (PF)  is  based  on  variational 
methods  and  the  theory  of  rearrangements  of  functions.  For  each  problem  it 
will  be  shown  that  there  is  a  critical  point  $  of  a  suitable  functional  and 
that  the  function  <t>  is  a  solution  of  a  weak  formulation  of  the  problem. 
Moreover,  ♦  can  be  taken  to  be  even  in  x,  nonnegative,  and  nonincreasing 
on  0  <  x  <  k  for  each  y.  That  is,  ♦is  symmetrized,  according  to  the 
following  definition. 


Definition  3.1.  Let  ♦  =  $(x,y)  be  continuous  on  [-k,k]  x  [0,1]  and  for 
each  y  let  y($,a,y)  denote  the  Lebesgue  measure  of  the  set 
(x  |  $(x,y)  >  a).  A  function  $  which  is  even  in  x,  nonnegative, 
nonincreasing  on  0  <  x  <  k  for  each  y  and  satisfies 

U(4>,a,y)  =  P(  |  <M  ,a,y) 

A 

is  called  the  symmetrization  of  If  ♦  =  ♦  we  call  it  symmetrized. 
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The  construction  of  $  and  the  properties  of  $  that  we  use  are  given 
in  an  appendix.  While  the  use  of  symmetrized  functions  is  not  needed  for  the 
existence  of  periodic  solutions,  it  plays  a  central  role  in  the  estimates  in 
section  4  which  we  use  to  show  that  a  solitary  wave  is  obtainable  as  the  limit 
of  periodic  waves  with  increasing  period.  A  symmetrized  function,  moreover, 
has  an  obvious  appeal  in  that  observed  permanent  waves  exhibit  such  a  form. 

We  note  that  there  are  solutions  of  (P)  of  any  period,  namely  the  zero 
solution  and  a  countably  infinite  collection  of  non-zero  x-independent 
solutions  obtainable  from  bifurcation  theory  (Crandall  and  Rabinowitz  1970, 
Turner  1973).  One  of  this  latter  class  of  solutions,  a  positive  one,  is  the 
conjugate  flow  studied  by  Benjamin  (1971).  We  shall  show  in  section  4  that 
the  solutions  we  obtain  here  are  x-dependent  provided  k  is  sufficiently 
large . 


The  analysis  will  be  carried  out  principally  in  the  Hilbert  space 
(R  *  (0,1))  defined  as  followst  let  be  the  space  of  C 

functions  which  have  support  whore  0  <  y  <  1  and  which  are  periodic  in  x 
with  period  2k.  Define 

lu»2  “  //  tlVuI2  +  h  (y)u2}dxdy  ,  (3.5) 

“k 


where 


*  {(x,y)  *  |x|  <  k,  0  <  y  <  1)  . 


The  Poincare  inequality. 


//  |Vu|  3xdy  >  w2//  u2dxdy  , 


(3.6) 


(3.7) 


together  with  the  inequality  hQ(y)  >  e  >  -w  ,  a  consequence  of  hypothesis 

All,  show  that  I  I  ^  is,  in  fact,  a  norm.  We  let  be  the  completion  of 

00 

in  the  norm  (3.5).  Then  is  a  Hilbert  space  and  the  corresponding 


0 
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inner  product  is 


F1 

i 

t 

*  i 

, 

; 

’  t 

i 


i 


(u,v)  =  //  {7u#Vv  +  h  (y)uv}dxdy 

*  a 

k 


(3.8) 


If  l  is  a  continuous  linear  functional  on  H^,  its  value  at  u  is  denoted 
by  < l , u> . 

We  will  also  have  occasion  to  use  the  spaces  ,  C-^'u,  L P,  and  W-*  'P  for 
various  domains  (cf.  Adams  1975,  Gilbarg  and  Trudinger  1977).  A  further 
consequence  of  the  Poincare  inequality  (3.7)  is  that  for  u  e  H^, 

lul.  <  C  //  |Vu|2dxdy  <  C  lul  .  <  C  lull  ,  (3.9) 

k  1  SL  2  W1/2(0.)  3  K 

k  k 


with  ,  C2,  and  C3  independent  of  k. 

In  what  follows  we  use  the  letter  C,  with  or  without  subscripts  or 
superscripts,  to  denote  a  positive  constant  which  may  be  different  in 
different  inequalities  and  which  may  depend  upon  pertinent  parameters,  but 
which  is  independent  of  k.  Moreover,  we  will  often  use  the  convention 

/g  *  //  g(x,y)dxdy  .  (3.10) 


Combining  (Adams  1975,  lemma  5.14)  with  (3.9)  one  finds  that  is 

compactly  embedded  in  2  <  p  <  "  and  that 

lul  <  Clul  ,  2  <  p  <  «  .  (3.11) 

LP<0)  k 

k 

If  u  e  H,  n  W2'*‘((i  )  and  v  6  H.  ,  then 
k  k  k' 

(u,v)^  *  /  (-Au  +  hQ(y)u)v  (3.12) 

the  boundary  terms  on  x  *  ±k  cancelling  because  of  periodicity.  We  are  thus 
led  to  define  a  weak  periodic  solution  of  (P)  as  a  pair  (X,$)>  <J>  e  H^, 
satisfying 

( '^ ' v ) vc  +  ^  =  *  /f(y ,4>)v  (3.13) 
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i 


i 


for  all  v  e  1^.  The  integrals  in  (3.13)  exist  by  virture  of  (3.11)  and 
assumptions  AI  and  All.  Weak  solutions  of  problem  (P)  will  be  obtained  from 
variational  principles  suited  to  the  separate  cases  (PC)  and  (PF). 

Consider  first  the  problem  (PC)  and  adopt  the  notation  (cf.  3.2),  that 
for  a  function  u  defined  on 


and  for  R  >  0, 


I(u)  -  J{|Vu|2  +  H(y,u)} , 

J(u)  *  / F(y,u),  \ 

S(R)  =  {u  e  H.  s  l(u)  -  R2}. 

k  ^ 


The  assumptions  on  f  and  h  guarantee  that  the  functionals  I  and  J  are 
defined  on  Hk-  One  calls  ♦  a  critical  point  of  J  on  S(R)  if  the 
derivative  of  J  in  directions  tangent  to  S(R)  is  zero;  that  is,  the 
derivative  J'(<fr)  is  parallel  to  !'(♦).  A  consequence  of  this  is  equation 


(3.13). 


•  ' 


Theorem  3.2.  Suppose  f  and  h  satisfy  condition  A.  Then  for  each 

2 

k  >  0  the  problem  (PCK3.3)  has  a  solution  (X^,^),  $k  e  Hk n  C  ^k^' 
satisfies 


1) 

J(4>k)  = 

sup  J( u ) , 
ues(R) 

2) 

O 

A 

•< 

and  >  0  in  fl, , 

k  k 

3) 

<t>  =  ■$ 

yk  *k 

(cf.  definition  3.1). 

• 

Proof .  We  fix  k  and  suppress  it  in  writing  X  and  $.  One  readily 
verifies  that  I  and  J  are  continuously  differentiable  and  that,  for  u, 
v  e  Hk, 
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<I'(u),v>  -  2/(7 u*Vv  +  h(y,u)v), 

<J'(u),v>  =  2/f(y,u)v. 

Moreover,  I'  is  strongly  monotone.  For,  using  the  Poincare  inequality 
(3.7),  we  see  that 

<I'(u)  -  l'(w),u-w>  =■  2/{|V(u-w)|2  +  (h(y,u)  -  h(y,w) ) (u-w)} 

>  2/{|7(u-w)|2  +  e(u-w)2} 

>  2(1  +  ~)  /|V(u-w)|2, 


•V  /V  2 

where  e  =  min(e,0).  Since  1  +  e/ir  >  0,  according  to  hypothesis  All,  it 
follows  from  (3.9)  that 

<I'(u)  -  I'(w),  u-w>  >  Clu-wl2.  (3.15) 

Then  from  the  standard  relationship 

I(u)  -  /g  <1* (Tu),u>dx 

and  (3.15)  one  has 


l(u)  >  C /’  TluI2dT 


(3.16) 


C  .  ,2 

2  ,UV 


Thus  all  functions  in  a  level  set  S(R)  of  I  are  uniformly  bounded  in  H^. 
Then  because  of  (3.11)  and  the  growth  conditions  required  of  f  in  hypothesis 
AI,  for  each  R  >  0,  J,  in  turn,  is  uniformly  bounded  on  S(R).  For  a 
fixed  R  >  0,  let 

sup  J(u)  *  b  . 
ues(R) 


Assumption  Al  implies  further  that  J(u)  =  0  if  and  only  if  u  *  0,  so 
b  >  0. 
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The  continuous  functions  which  are  piecewise  linear  in  x  and  y  and 

have  nonzero  partial  derivatives  with  respect  to  x  a.e.,  form  a  dense  linear 

subspace  of  (cf.  the  appendix) ,  so  there  is  a  sequence  of  such  piecewise- 

linear  functions  Uj  e  S(R),  j  «*  1,2,...,  for  which  J(u^)  ♦  b  as  j  ♦  “. 

Letting  Uj  be  the  symmetrization  of  Uj,  we  see  from  relation  (Ap2)  of  the 

appendix  that  J(Uj)  ■  J(Uj)  for  j  -  1,2,...,  while  from  (Ap2)  and  (Ap3) 

l(u. )  <  I(u  )  ■  R2.  Since  lu  I  is  uniformly  bounded,  according  to 

]  J  J  K 

(3.16),  u^  has  a  subsequence  converging  weakly  to  a  function  $  e  H^.  We 

will  show  $  is  a  solution  of  (PC). 

As  was  noted  in  conjuction  with  inequality  (3.11),  Hk  is  compactly 

embedded  in  Ln+^(0.),  where  n  occurs  in  hypothesis  All,  so  that  the 

subsequence  converges  strongly  to  $  in  Ln+ 1 ( )  •  Moreover,  a  further 

subsequence  of  {(L)  may  be  assumed  to  converge  pointwise  a.e.  to  $  on 

According  to  (Krasnoselskii  1963,  p.  30)  the  map  u  ♦  F(y,u)  is 

continuous  from  Ln+1(ft^)  into  L^fl^).  Thus  J($)  *  b,  and  as  b  >  0, 

HO.  The  functional  l  is  the  sum  of  the  norm  in  H^,  which  is  lower 

semi- continuous,  and  the  term  JH^ (y,u)  which,  like  J(u),  is  continuous. 

2  2 

Hence  I($)  <  R  .  Suppose  that  I($)  <  R  .  Since  I*  is  strongly  monotone, 

the  function  I(t<fr)  is  increasing  in  t  and  approaches  +  *  as  t  ♦  +  “. 

Since  I(t<fr)  is  also  continuous  in  t,  there  is  a  t  >  1  for  which 
2 

I(t$)  =*  R  .  Hypothesis  AZ  implies  that  F(y,u)  is  strictly  increasing  in  the 

variable  u,  so  J(t$)  >  b.  This  contradicts  the  characterization  of  b  and 
2 

hence  I($)  =  R  .  Since  <!'(♦)»♦>  >  0  we  can  use  a  lemma  of  Lyusternik 
(Vainberg  1964,  p.  96)  to  conclude  that  J’($)  is  a  multiple  of  !'($)• 

Since  J'($)  ?  0  we  can  write  this  as 

!'(♦)  -  VJMH, 


and  this  is  tantamount  to  (3.13). 
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It  follows  from  hypothesis  AI  that  <J‘ ($),♦>  >  0  and  hence  X  >  0.  As 

regards  $,  since  it  is  a  weak  solution  of  an  elliptic  equation,  it  follows 

from  the  elliptic  theory  (Agmon,  Douglis  and  Nirenberg  1959)  and 

inequality  (3.11)  that  is  in  W2'^^)  for  any  p  <  •.  The  Sobolev 

1,o 

embedding  theorem  shows  that  ^  e  C  (fl^)  and  the  Schauder  theory  (Gilbarg 
and  Trudinger  1977)  yields  ♦  e  C2,°(fl^).  Being  the  pointwise  limit  of 
symmetrized  functions,  the  continuous  function  $  inherits  the  property  of 

A 

being  symmetrized;  i.e.  $  *  $.  To  see  that  $  is  strictly  positive  on  £1^ 
we  use  hypothesis  All  to  write  h(y,$(x,y))  =  h(y,$(x,y) )$(x,y) ,  and  then 
decompose  h  as  h  (x,y)  -  h  (x,y)  with  h  being  nonnegative  functions 
of  x  and  y.  Since  ♦  >  0,  -A$  +  h+$  =  h  <p  +  Xf(y,$),  a  nonnegative 
function  which  is  not  identically  zero.  The  strong  maximum  principle  (Gilbarg 
and  Trudinger  1977)  implies  that  $  >  0  in  £1^,  completing  the  proof. 


I 


We  next  turn  to  the  problem 


(PF). 


For  use  here  and  in  the  sequel. 


denote  the  smallest  eigenvalue  of  the  linear  problem, 

-Au  +  hQ(y)u  «  XfQ(y)u, 


I  (3.17) 

u  e  Hk, 

by  p.  One  obtains  the  eigenfunctions  of  (3.17)  by  separation  of  variables 
and  it  is  easily  seen  that  |i  has  an  associated  eigenfunction  5(y), 
independent  of  x.  It  is  the  eigenfunction  associated  with  the  lowest 
eigenvalue  of 


"Cyy  +  ho(yK  =  Xfo(y)5 


(3.18) 


5(0)  =  5(i)  =o  ) 

and  is  known  to  be  positive  (Ince  1927,  p.  235).  The  existence  and  positivity 
of  5  also  follow  easily  from  the  type  of  arguments  just  completed  for  the 
problem  (PC). 
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It  will  be  convenient  to  normalize  5  so  that 


f'o  +  V2)  =  1-  (3>19) 

For  fixed  X  define  a  functional  M  on  by 

M(u)  -  lul3  +  Ja  {h^ (y,u)  -  XF(y,u)}.  (3.20) 

k 

We  use  lM’(u)l  to  denote  the  norm  of  the  derivative  M'(u)  as  a  functional 
on  H^.  One  readily  verifies  that  a  critical  point  of  M,  i.e.  a  point  $ 
for  which  <M'($),v>  **  0  for  all  v  e  H^,  is  a  weak  solution  of  problem  (PF) 
(cf.  3.13).  To  show  there  is  a  critical  point  we  use  a  slight  modification  of 
a  result  of  Ambrosetti  and  Rabinowitz  (1973,  theorem  2.1). 


Proposition  3.3.  Let  k  >  0  be  fixed  and  let  M  be  a  continuously 
differentiable  functional  on  satisfying  the  following  conditions. 

1)  There  are  constants  a  >  0  and  r  >  0  such  that  M(u)  >  a  for 

lul  ■  r. 
k 

2)  M(0 )  “0  and  there  is  a  function  w  with  Iwl,  >  r  such  that 

k 

M(w)  <  0. 

3)  For  each  0  >  0  if  is  a  sequence  satisfying 

B  <  M(u^)  <  0  \  for  all  i,  and  lM'(u^)l  ♦  0,  as  i  +  •»,  then  a 

r  ,  «  m 

subsequence  tu^  of  ‘ui^i=i  converges  strongly  in  as 

m 

m  ♦  *. 


Let 


and  define 


r  -  tY  e  c([o,i],Hk>  |  y (o )  =  o,  y(D  »  w> 


b  =  b  (M)  =  max  M(u) 
1  ueY([0,D) 


b  *  inf  b  . 

Yer  ' 


(3.21) 


Suppose  Y  ,  n  =  1,2,3,...,  is  a  sequence  of  paths  in  T  such  that 


•  ' 


•  • 
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Vb*"- 


,  ,  m 

Then  there  is  a  subsequence  ln^)  ^  of  the  positive  integers  and  elements 

u  e  Y  ([0,1])  such  that  uR  converges  strongly  to  u  as  j  ♦  “, 
nj  Uj  j 


M(u)  =  b,  and  M'(u)  =  0. 


Proof.  For  a  real  number  z  introduce 


A  =  {u  e  H.  |  M(u)  <  z} 
z  k 


Kz  =  {u  e  |  M(u)  *  z  and  M*  (u)  -  0) . 

For  j  =  1,2,3,...  let 

N.  =  {u  e  A  \A  |  IN'  (u)  I  <  i). 

3  ”"7 

Then  for  each  j,  Nj  contains  a  neighborhood  of  K^.  From  (Ambrosetti  and 
Rabinowitz  1973,  lemma  1.3)  it  follows  that  for  each  j  there  is  a 
deformation  map  n^(x),  taking  (t,x)  e  [0,1]  *  continuously  into 

and  constants  ej  >  £j  >  0  such  that 


>  0 

such 

that 

1) 

no(u) 

a  u 

for  u  6  H^, 

2) 

nj(u) 

■  U 

for  u  e  A  , 

b-‘J 

3) 

"j’  '  Ve.- 

For  each  j  choose  a  positive  integer  nij  so  that  nij  <  £j* 

Y  ([0,1])  n  N.  =  0,  then  since  the  composition  h^oY  is  a  member  of  r 
mj  J  1  “j 

and  M(u)  <  b-£j  on  the  image  of  this  composition,  we  have  a  contradiction  to 
the  characterization  (3.21).  Hence  there  is  a  function 

u  e  Y  ([0,1])  n  N,  for  each  j.  The  compactness  property  (3)  guarantees 
m  m .  j 

3  3 


I 


that  a  subsequence  {u  }  of  tu  }  converges  strongly  as  j  ♦  •  to  a 

n .  m . 

j  3 

function  u,  from  which  it  follows  that  M(u)  *  b,  and  M'(u)  =  0.  Note 
that  the  hypotheses  of  the  proposition  imply  that  b  >  a  >  0  and  thus  u  f  0. 

We  can  now  show  that  problem  (PF)  has  a  symmetrized  solution. 

Theorem  3.4.  Let  f  and  h  satisfy  condition  A  with  a  and  0  as 

defined  therein.  Let  y  be  the  lowest  eigenvalue  of  problem  (3.17).  If 

a  <  X  <  vi  then  for  each  k  >  0  the  problem  (PF)(3.4)  has  a  solution 
2 

♦k  6  n  C  (8)  with  the  following  properties. 

1)  M(4>  )  =  inf  b  (M)  (cf.  3.20  -  3.21), 

*  YeT  T 

2)  $k  >  0  in  8, 

A 

3)  <|>k  =  $k  (c.f .  definition  3.1). 

Proof .  Fix  k  >  0.  Since  y  is  the  smallest  eigenvalue  of  (3.17)  it 
follows  that 

lul^  >  y/f0(y)u2  (3.22) 

for  any  u  e  H^.  Hypothesis  AI  implies  that  for  any  6  >  0  there  is  a 
constant  >  0  so  that 

If t (y»t) |  <  6|t|  +  Cj |tn| , 

from  which  it  follows  that 

/^(y#  u)  <  /  C«5u2  +  c£|u|n+1j. 

Upon  invoking  the  embedding  inequality  (3.11)  we  find  that 

/f^ (y,u)  <  <C6  +  c£  lulk-1)«ul2, 

so  that  the  functional  / F,  is  o(lul2)  as  lul,  ♦  0.  One  obtains  a 

'1  k  k 

corresponding  estimate  with  F^  replaced  by  H^.  These  estimates,  with  the 
inequality  (3.22),  lead  to 


M(u)  >  (1  -  ^lul^  + 

as  lul^  *  °*  So  f°r  ®ach  A  <  m,  M  satisfies  hypothesis  (1)  of  proposition 
3.3  for  any  sufficiently  small  positive  r. 

Plainly  M(0)  =  0.  To  find  a  function  w  such  as  occurs  in  hypothesis 
(2)  of  proposition  3.3,  first  observe  that  from  AI  and  All, 

Af^ (y»t)  -  H^yjt)  >  (Aa-a)tm+1/(m+1). 

Now  if  u  is  any  positive,  symmetrized  function  in  and  s  >  0,  then 

r».  ,  .  ..  ,  .  (Ao-a)sm+1  r  m+1 

J  (Af^ (y,su)  -  H1<y,su)}  >  - — -  Ju  , 


and  so  from  (3.20), 


m+1 


...  2..  .2  y  f.  ,  .  2.  (Aa-O)s  r  m+1 

M(su)  <  a  (lu*k  “  AJ f  (y )u  ) - — -  Ju 

Then,  since  m  >  1,  for  s  sufficiently  large,  say  s  =  Sq,  M(squ)  <  0  and 

A 

ls0ulk  >  r.  We  take  w  =  sQu  and  have  w  =  w.  That  M  is  C1  and 

satisfies  hypothesis  (3),  often  called  the  Palais-Smale  condition,  is  shown  as 

in  (Aubrosetti  and  Rabinowitz  1973,  lemma  3.6).  Proposition  3.3  is  now 

applicable  and  the  existence  of  a  solution  of  problem  (PF)  satisfying  (1)  can 

be  concluded.  However,  establishing  a  solution  satisfying  (1),  (2),  and  (3) 

requires  further  argument.  To  that  end  it  will  be  verified  that  there  is  a 

sequence  of  paths  Y  e  F,  n=1,2,3,*»»,  with  images  consisting  solely  of 

n 

symmetrized  functions  for  which  b^  converges  to  the  critical  value  b. 

1  n 

Proposition  3.3  will  then  yield  a  critical  point  satisfying  properties  (1), 
(2),  and  (3). 

We  begin  by  showing  that  the  functional  M  is  Lipschitz  continuous  on 
bounded  subsets  of  H^.  Let 

G1 (y,t)  =  AF(y,t)  -  ( y, t)  (3.23) 


and 


Then  from  (3.20), 


N(u)  =  Jg  (y,u). 


-24- 


i 


< 

1 

\ 


I 


|M(u)  -  M(v)  |  <  |«u«k  "  ,vl^ 

+  |/q  <N'(x(u-v)  +  v,  u-v>  dx  | 

-  (Iulk  +  lvlk)  |lulk  -  lv«kl 

+  |/q  /fl  [g^y^lu-v)  +  v]  (u-v)dx  | 
k 

<  (lul.  +  I  vi .  +  c  fl  ■  g, (y >x  (u-v)  +  v)l  dx  )  I  u-vl  , 

k  k  0  1  L  <V 

where  gj(y,t)  “  a^Gjty.t).  The  growth  conditions  in  AI  and  All  imply  that 

the  expression  in  parentheses  is  bounded  on  any  bounded  subset  of  H^. 

Let  Y  be  a  path  in  T  and  suppose  lY(t)lk  <  s  for  0  <  t  <  1.  Let 

B  denote  the  ball  of  radius  s  centered  at  0  in  Ht  and  let  C  denote 
s  K 

the  Lipschitz  constant  for  M,  associated  with  the  ball  •  Given  e  >  0 
(e  <  4Cs)  choose  a  partition  0  *  tg  <  tj  <• • • <  tn  =  1  so  that 
ly(ti+1)  -  T(ti)»k  <  e/4C  for  i  ■  0,1,...n-1.  Let  uQ  «  0,  un  *  w,  and  for 
i  -  1,2,»»»,n  let  uA  be  a  piecewise  linear  function  of  x  and  y  having 
3  u.  r1  0  a.e.,  such  that  I  u.  -  Y(t.  )l,  <  e/4C  (cf.  the  appendix).  Define 

XX  1  1  X 

a  path  Y  by  linearly  interpolating  the  functions  u^j 


^  1  {+  1 
Y  (t)  -  (1  -  (nt-i))^  +  {nt-i)ui+1,  J  <  t  < 


i  »  0,1,2, ••• ,n-1. 


(3.24) 


Since  e  <  4Cs  each  u^  is  in  Bjg  and  hence  the  path  Y  lies  in  B2g. 

Using  the  Lipschitz  property  of  M  one  checks  that  b^  <  b +  e/2.  Suppose 

Y  T 

#  u  ■  —  <  Q  for  i  »  0,1, Then  for  each  t  €  [0,1],  Y(t)  is  a 

L  Bk) 

convex  combination  of  some  pair  u^,  ui+1  and  l*ence  lY(t)>  m  <  Q. 

L 

It  will  now  be  convenient  to  let 

G(y,t)  =  *F(y,t)  -  H(y,t) 


and  express  M  as 
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M(u)  *  /{|7u|2  -  G(y,u)). 

Since  G(y,t)  is  uniformly  continuous  for  0  <  t  <  Q,  there  is  an  h  >  0  so 

that  if  |t-t|  <  n,  then  |G(y,t)  -  G(y,t)|  <  e/4k,  where  2k  is  the 

period.  We  suppose  that  lu^  -  u^+^l  <  h  for  i  *  0,1, ...n-1.  If  this  is 

L 

not  the  case,  the  t  interval  may  be  partitioned  more  finely  so  that  at 

~  m 

adjacent  points  of  the  partition  the  values  of  Y(t)  differ  in  L  by  at 

a# 

most  h.  Assuming  that  Y  is  defined  by  (3.24)  for  a  possibly  finer 

A 

partition,  let  Y  be  the  path  obtained  by  replacing  each  function  u^  in 

A  A  A 

(3.24)  by  its  symmetrization  u^.  A  convex  combination  of  u^  and  u^+1 

A 

will  be  symmetrized  and  so  Y  consists  entirely  of  symmetrized  functions. 

A  A 

Note  that  w,  the  end  of  the  original  path,  satisfies  w  *  w.  The  path  Y 
is  a  substitute  for  the  set  one  would  get  by  symmetrizing  each  point  Y(t)  on 
the  path.  Symmetrizing  each  point  on  a  path  does  not  obviously  produce  a  path. 

*  A.U 

Now  consider  M  on  the  path  Y.  On  the  i  interval  of  the  partition 

A  A  A 

used  in  (3.24),  Y(t)  may  be  expressed  as  (1  -  Tju^  +  Tui+i  where 
T  =  nt  -  i.  By  convexity, 

/  |V((  1-T)U±  +  TUi+1)|2  <  (1-T)/|7ui|2  +  t/  l?Ui+1|2. 

Next,  using  inequality  (Apl)  of  the  appendix  it  is  concluded  that 


and  so 


,ui  •  Vi'  -  <  ,ui  "  ui+i'T-  <  n* 

L  L 


A  A  A  A 

/ I G(y , ( 1-T )uA  +  T«i+1)  -  (1-x)G(y,ui)  -  TG(y,ui+1)| 

A  A  A 

<  (1-T)/|G(y,(1-T)ui  +  Tui+1>  -GCy,^)! 

A  A  A 

+  T / |G(y,(1-T)Ui  +  Tui+1)  -  G(y*ui+t)l 

<  /  — 

1  4k 


because  of  the  uniform  continuity  of  G.  Finally,  using  relations  (Ap2)  and 

,  A 

(Ap3 )  of  the  appendix  we  estimate  M  on  the  1th  segment  of  Y  by 
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M( ( 1-T )ui  +  TUi+1)  <  ( 1-T )/ | Vu±|  +  T/|7ui+1l 

-  (I-Tj/Gty,^)  -  t/G(y,ui+1)  +- 

-  ( 1-T)M(ui>  +  TM(ui+1)  +  | 

<  (I-tJmCu^  +  TM(ui+1)  +  | 

<  b  +  | 

7  2 

<  by  +  e. 

A 

As  the  inequality  holds  on  each  segment  of  Y  , 

b»  <  b  ♦  e. 

Y  T 

If,  for  any  positive  integer  n,  Y  is  a  path  for  which  b  <  b  +  — , 

n  T  n 

n 

according  to  what  we  have  just  shown,  there  is  a  path  Yn  consisting  of 

2 

symmetrized  functions  for  which  b  <  b  +  —  •  Then  from  proposition  3*3  there 

^n 

m  A 

is  a  sequence  {v.},  .  with  v.  **  v , ,  for  all  j,  converging  in  H,  to  a 

j  j  3  K 

nonzero  critical  point  of  M.  It  now  follows  exactly  as  in  the  proof  of 

_  A 

theorem  3.2  that  is  in  C2'  (fl^),  that  and  that  >  0  in 


Remark.  The  restriction  X  <  y  in  theorem  3.4  stems  from  our  ultimate 
interest  in  waves  of  supercritical  speed  which,  as  will  be  seen,  persist  in 
the  limit  of  increasing  periods.  While  one  can  show  the  existence  of 
nontrivial  periodic  solutions  of  problem  (P)  having  X  >  y,  using 
bifurcation  theory,  they  are  not  useful  for  our  present  study.  Assuming 

X  <  ii,  the  further  restriction  o  <  X  is  essential  as  one  can  see  from 

3  3  2 

the  example  +  0$  =  X($  +  a$  )  in  which  y  =  w  .  Suppose  there  were  a 

nontrivial  solution  (X  ,$ )  with  X  <  a  ^3.  Taking  the  inner  product  of  each 
side  of  the  equation  with  $  one  finds  that  / |7$|2  <  X/$2  and  hence  from 
(3.7),  X  >  n  ,  which  is  excluded  by  the  restriction  X  <  y. 


One  can  see  from  the  proof  of  theorem  3.4  that  it  is  the  behavior  of  the 
whole  term  Xf^  -  h^  that  is  important  in  problem  (PF).  Thus,  conditions  AI 
and  All  could,  in  this  problem,  be  replaced  by  an  analogous  condition  on 
Xf  -  h,  so  allowing  for  wider  applicability  of  the  theory. 


4.  BOUNDS  FOR  PERIODIC  SOLUTIONS 


%  * 

The  subsequent  analysis  depends  on  bounds  on  the  solutions  (Ak,$k) 

(3.1)  whose  existence  was  established  in  section  3.  In  the  problem  (PCH3.3) 
the  size  of  l$l  was  specified  while  in  the  problem  (PF)(3.4)  A  was 

X  K  £ 

specified  in  the  range  a-1<J  <  A^  <  y.  It  will  be  important  to  have  estimates 
independent  of  k  for  both  A^  and  $k  the  two  problems.  Such  bounds 
will  be  obtained  in  the  present  section,  using  some  additional  hypotheses 

'  •  ' 

which  are  not  unduly  restrictive  for  the  physical  applications. 

It  will  first  be  shown  that  A^  <  y  in  the  problem  (PC).  When  it  is 
convenient  the  subscript  k  will  be  suppressed. 


Lemma  4.1.  Suppose  that  f  and  h  satisfy  condition  A  of  section  3, 
with  the  additional  restrictions  m  <  5  and  2a  *0  <  y  on  the  parameters. 

It  follows  that  there  are  constants  n  *  MR)  >  0  and  kQ  -  kQ(R)  such  that 
if  *-8  th®  solution  of  (PC)  obtained  in  theorem  3.2,  with 

I(<j»k)  -  R2,  then 

Ak  <  (1-n)V,  (4.1) 

for  k  >  kg(R). 

Proof .  We  suppress  the  variable  y,  writing  f  (♦)  for  fj(y,$),  and 
so  on.  In  the  weak  form  (3.13)  of  the  problem  (3.1),  let  v  *  ♦  *  $k  and  use 
the  notation  (3.14)  to  express  A  -  Ak  as 

/{ Wi2  +  hQ^2  +  »,(+)*} 

/{ fQ$2  +  «,<♦)♦) 

R2  +  /{hn (♦  )«|»  -  H t (4 ) } 

=  J<*)  ♦  J  <f 1 <♦)♦  -  F,($)}  * 
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Since  J(u)  <  J($)  for  an  arbitrary  u  e  with  I(u)  “  R2,  then 


R2  +  /{h1<^)*  -  H1(«fr )} 

X  <  J(u)  +  J {£ t(4> >♦  -  F^*)}* 


(4.2) 


Take  u(x,y)  to  be  the  element  of  defined  on  -k  <  x  <  k  by 

u(x,y)  -  A(x)£(y),  where  5  is  the  eigenfunction  defined  in  (3.18).  We  show 

2 

that  A  can  be  chosen  so  that  I(u)  «  R  and  so  that,  for  k  >  kQ,  the 

right-hand  member  of  (4.2)  is  at  most  (1  -  n),  where  kQ  and  n  are 

positive  numbers  depending  on  R. 

2 

The  normalization  I(u)  «*  R  entails  having 

R2  -  n2  /  J  A2  +  J*k  A2  +  fQ  H1  (u) ,  (4.3) 

k 

where,  for  0  <  p  <  00 , 

nP  *  n  ?Pdy*  (4*4> 

r  1  2  2 

Here  we  have  used  the  normalization  J  0^y  +  h^C  }  *  1  given  earlier  (cf. 
3.19).  The  normalization  also  implies  J^f^C2  “  W  \  and  thus 

J(u)  -  U_1{/^K  A2  +  w/n  F1 (u) } .  (4.5) 


Combining  (4.2),  (4.3),  and  (4.5)  one  sees  that  the  inequality  (4.1)  will  hold 
provided 


n2  J-k  Ax  +  n  f-k  ft2  <  /{  (1-n)HF|(u)  -  Ht(u)> 

+  /{ ( 1-n  )m  (f  t  -  f^  ($ ) )  -  (h^*)*  -  H1  (♦ ) )} . 


(4.6) 


From  condition  AIII  (in  this  instance  0  <  1  suffices)  it  can  be  seen  that 
( i-n )u (f , (♦  )♦  -  f^))  -  (hn<^ )♦  -  h  1  ($ ) ) 

>  (i-n)uf  ,(+)*  -  h^*)* 

-  s ( ( i-n )Mf  1  <♦  )♦  -  h^^)^) 

-  (i-0)((i-n)yf  (+)+  -  h^*)*) 


‘  • 
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>  ( i-e ){ < i-n )wa  -  o}*»+1. 

Hence  it  suffices  to  find  an  0  with  0  <  n  <  1/2  so  that, 

rk  .  2  .  _rk  .2 


n2  ^-k  x  *  n^-k  A  *  J  t(1’n)w  F!(u)  "  Vu)1 


(4.7) 

-JC  1  I 

for  some  u  satisfying  (4.3). 

To  continue,  let  A(x)  -  R8g(\>x)  where  g(t)  -  e-2^  and  where  8  and 
v  are  constants  to  be  determined.  A  short  computation  brings  (4.3)  to  the 
form, 

,  ft2  am+1 

1  -  4n2V0  c2(Vk)  +  c2(Vk)  +  cm+1(v>0,  (4.8) 

where  c  (s)  ■  JB  gm(t)dt  and  where  5  -  6(k)  is  a  constant  satisfying 
m  “8 

-2J'n  <  (m+1)6  <  20  n  . ,  (cf.  condition  All  in  section  3).  Using 

m+ 1  ntT  I 

conditions  AI  and  All  one  sees  that  (4.7)  will  hold  provided 

2 


4n2v02c2(vk)  +  n  c2(vk)  <  2n. 


m-!  “+1 

<4-5’ 


where  Y  ”  (l-h)pa  -  0  >  0.  If  we  now  let  v  «  8 


2+e 


and  require  n  to 


satisfy  the  further  condition  h  <  n2&4+2E  then  (4.8)  and  (4.9)  will  follow 


from 


and 


1  ■  4n^4+€c2  +  ®  e®2  +  1  C° 


m+1 ' 


5n34+€c  4  2n 

2  2  m+1 


ra-1»m-1-e 
YR  8  c 


m+1 


m+1 


(4.10) 


(4.11) 


2+e, 


respectively.  Note  that  cp(p  ■  2, m+1)  now  has  8  k  as  its  argument  and 
further,  that  as  the  argument  increases  to  infinity,  c^  increases 
monotonically  to  1/p.  We  will  henceforth  assume  that  k  >  k(8)  * 

(^  ln2)8  2-£  so  that  1/2p  <  cp  <  1/p. 

To  satisfy  (4.11)  it  will  suffice  to  have  8  <  8^(R)  where 


S^R)  *  ( 


2nm+  1^R 
5n2(m+1 ) 


m-1 


1 5-m+  2E 


(4.12) 
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Here  we've  used  the  bound  cm+j/c2  *  1/(m+1).  To  satisfy  (4.10)#  begin  by 

"1/eo 

choosing  an  with  0  <  Eg  <  m-1  and  letting  8g  *  4  ,  so  that 

0O  °c2  >  1.  Our  choice  for  8  is 


0(R) 


inf  min 
0<e<eo 


1 


(m+1)' 


80n  R 
m+1 


Si) 


m-  1-e 


,  eo,01(R) 


He  further  choose 

kQ(R)  “  ln2)(0(R))  ° 

and 

4+2£ 

0{R)  *  min{-,  n2(8(R))  >• 

-1/eo 

Since  0Q  *  4  <  1,  8  <  0Q  <  1  and  so  for  all  e  in  [0,eq), 

4+  v 

n  <  njt.  ,  as  required.  Assume  k  >  kQ *  Then  for  e  in  [0,Eg], 
k  >  kg  >  k(0),  from  above,  so  1 /2p  <  cp  <  1/p. 

Since  0(R)  <  01<R),  (4.11)  is  satisfied.  By  the  choice  of  0(R)  we 

also  guarantee  that  the  first  and  third  terms  in  the  right-hand  side  of 
(4.10)  are  at  most  1/4  for  any  £  in  [0,Eg].  The  remaining  term,  0  E(R)c2 
approaches  c2  <  1/2  as  e  ♦  0  and,  since  0(R)  <  0Q,  exceeds  1  for 
£  -  Eg.  By  continuity  there  is  an  e  in  [0,Eg]  for  which  (4.10)  is 
satisfied  and  we  make  it  our  choice.  The  trial  function  u  is  now  fixed  and 
the  inequality  (4.1)  follows  accordingly. 


Remark.  As  R  approaches  0,  one  has 

«  (m-1)(4+2e0) 

n(R)  =  0(R  ),  where  oi  =*  - - - 


Lemma  4.2.  Suppose  that  f  and  h  satisfy  condition  A  of  section  1.  Then 
there  are  constants  Cj  and  C2  such  that  for  any  solution  ( X , <t> )  of  (PC) 
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- — — - ' : ’-1 

(4.14) 

Proof .  As  in  the  proof  of  lemma  4.1 

1*1*  +  /  »»,<♦>♦ 

X  - - - - . 

/{*„♦  +  «,(♦)♦) 

2  (  2  ® 
Since  1*1*  >  y/f0^  by  (3.22), 

1  +  i<M“2/h1(^)<fr 

x  >  _ *  1 _ 

-1  -2t  ' 

V  +  l*«k  /*,(♦>♦ 

and  the  desired  estimate  follows  from  the  growth  assumptions  on  f 1  and  h1 
together  with  inequality  (3.11). 

We  now  have  upper  and  lower  bounds  for  X  occurring  in  the  solution 
(*,♦)  of  problem  (PC)  given  in  theorem  3.2.  In  the  proof  of  that  theorem  the 

estimate  l$l.  <  CR  (cf.  3.16)  was  also  derived.  In  the  problem  (PF)  X  is 

-1  * 
restricted  to  the  interval  (o  <x,w)  at  the  outset,  in  theorem  3.4.  For  this 

latter  problem  we  next  obtain  a  bound  for  11$!^  in  terms  of  X. 


with  *  e  Hk  and  l($) 


R  >  0  one  has  X  >  0  and,  moreover, 

1 

1  -  0* C  R 

x  >  - !srr)* 

1  +  ac2R 


Lemma  4.3.  Suppose  f  and  h  satisfy  condition  A  and  further  that 
a  ^  <  X  <  ii  and  1  <  m  <  5.  Then  the  solution  $  of  problem  (PF)  obtained 
in  theorem  3.4  satisfies 


<  C 


(V-X) 


5-m 
4 (m- 1 ) 


(Xa  -  a) 


(4.15) 


where  C  depends  on  3  and  u . 
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Proof.  He  will  use  the  variational  characterization  of  $  in  theorem 


3.4  to  obtain  (4.15).  Prom  proposition  3.3,  for  any  path  Y  €  I* , 

**'k  +  "  XF<y»*»  =  b  <  b^.  (4.16) 

On  the  other  hand  the  weak  form  (3.13)  of  the  basic  equation,  with  v  =  $, 
yields 

■  ♦I*  +  *  XV2  “  -  0*  (4.17) 

Combining  (4.16)  and  (4.17),  and  using  assumption  AIII,  one  arrives  at 

(1-8)(l$l2  -  X/f0(y)*  )2<  by.  (4.18) 

How,  by  (3.22), 

(1-8)<*l2(1  -  £)  <  bY,  (4.19) 

and  this  holds  for  any  Y  €  T . 

We  now  construct  a  particular  path  Yq  as  follows.  Let  u  be  the 
function  in  defined  by, 

u(x,y)  -  exp(-(y-X)1/2|x|)C(y), 

for  (x,y)  €  [-k,k]  x  [0,1],  where,  again,  5  denotes  the  eigenfunction 
associated  with  the  lowest  eigenvalue  of  (3.18),  normalized  as  in  (3.19). 

From  the  proof  of  theorem  3.4  there  is  inferred  the  existence  of  a  tQ  >  0 
such  that,  M(tgU)  =  0  where  M  is  the  functional  defined  in  (3.20).  Let 
YQ  be  the  path  defined  on  0  <  t  <  tQ  by  YQ(t)  *  tu.  A  simple  calculation 
shows  that 


M( tu)  < 


tW2 


ill  f052  + 


.  ,„..m+1  (l  -  exp[-2(m+1  )k(y-X)  ^2] ) 

+  (a-Xa)t  Vi 

(m+1 )  (y-X) 


(4.20) 


where  n^  is  defined  in  (4.4).  By  definition,  the  maximum  of  M(tu),  for 

t  in  the  range  0  <  t  <  tQ  is  b  ,  and  another  calculation,  using  (4.20), 

Yq 
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shows  that 


\  ‘ c' 


(u-X) 


m+3 

2(m-1) 


(Xa-o) 


2 

B-1 


(4.21) 


The  combination  of  (4.21)  and  (4.19)  with  Y  *  Yg  produces  the  desired 
inequality  (4.15). 


As  already  noted,  the  smoothness  of  a  solution  $  follows  from  standard 
regularity  theory.  We  repeat  some  of  it  here  to  exhibit  estimates  with 
constants  independent  of  the  period  2k.  Zn  particular  lemma  4.9  below  will 
be  used  in  showing  the  exponential  decay  of  solutions  in  lemma  4.10. 

CO 

Let  n  -  n(x)  be  an  even  C  function  supported  on  -3/4  <  x  <  3/4, 
taking  values  in  [0,1],  and  so  chosen  that  n(x)  +  0<1-x)  ■  1  for 
0  <  x  <  1.  This  last  property  is  equivalent  to  having 

n(x)  -  1/2  odd  with  respect  to  x  *  1/2  on  the  interval  0  <  x  <  1.  We  omit 
a  detailed  construction  of  such  a  function.  Given  n  we  let  n  (x)  -  n(x-s) 


for  any  real  s  and  note  that  the  collection  n j,  j  an  integer,  forms  a 

partition  of  unity  on  the  line.  Let  C  *  C(x)  stand  for  an  arbitrary  cutoff 

00 

function,  which  is  to  say,  C  is  a  C  function  with  range  in  [0,1]  having 

compact  support.  We  use  the  estimate, 

lul  <  C( lAul  +  lul  +  Bui  ^  ),  (4.22) 

W *,P(V)  LP(D)  LP(D)  L  OP) 


(Agmon,  Douglis,  Nirenberg  1959,  thm.  15. 2)  valid  for  a  domain  D  with  smooth 
boundary  3D.  The  constant  C  in  (4. 22)  will,  in  general,  depend  on  the 
domain.  However,  as  (4.22)  will  be  used  only  for  x-translates  of  certain 
fixed  domains,  that  dependence  will  not  enter.  Let 

S  .  “  { (x,y)  e  n  |  |x-a|  <  b),  (4.23) 

a,b 
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and  suppose  C  has  its  support  in  S 


Then  from  (4.22),  if 


a,b‘ 


u  e  w*;*  n  Hk, 


ICul  <  c(lA(Cu)J 

W  'P(Sa  )  ) 

a ,  b  a,  b 


+  I  Cull 


L  (Sa  h} 
a,  b 


),  (4.24) 


since  Cu  vanishes  on  3s  ,.  The  domain  s,  v  has  corners,  but  since  the 

a,b  a,b 

support  of  Cu  is  at  a  positive  distance  from  the  corners  the  estimate  can  be 
seen  to  hold.  Now,  since  A(Cu)  =  CAu  +  2Vc»Vu  +  uAC  it  follows  easily  from 


(4.24)  that 


ICul 


W2,P(S  > 
a,b 


<  C'(llAul 


LP<s.,b> 


«ul  .  ). 

W1,P(S  .) 
a,b 


(4.25) 


Lemma  4.4.  Suppose  f  and  h  satisfy  condition  A  and  let  (X,$)  be  a 
solution  of  problem  (P)(3.1)  having  0  <  X  <  y.  Then  for  each  integer  j, 

.  <  C(l*l  .  +  l*ln  ),  (4.26) 

o  <«,,,)  «  -2(Sj,3)  w,,2isj>3) 


where  C  depends  on  y  and  the  parameters  entering  condition  A. 


Proof .  The  growth  assumptions  in  AI  and  All  yield  an  estimate 

I  A*  |  <  C'(  I<t>|  +  |<|>|n),  (4.27) 

valid  at  each  point  (x,y),  with  C'  depending  on  y  and  the  parameters 

from  AI  and  All.  If  we  use  C  =  r'j_^  +  in  (4.25),  there  obtains 

a  W2,p  estimate  of  $  on  Sj^  in  terms  of  1^  estimates  of  A<{>  and  <}> 

on  ® j ,  2*  Estimates  on  ®j,2  an<^  3  are  similarly  related,  and  in  each 

estimate  the  constant  will  involve  the  fixed  function  n.  Combining  these 

estimates  with  the  Sobolev  embedding,  theorem  for  some  p  >  2,  one  has 

(letting  S^  .  =  S^ ) , 

J  /  1  J 
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»<M  .  <  C.I*I  ,  n  <  C_(dA4>*  +  »♦'  .  _  ) 

c 1  (s  1 )  1  w2,p(s1)  2  lp(s2)  W1,P(S2) 

<  C_(«*l  +  *♦!"  +  »♦«  2  ,  ) 

LP(S2)  lp(s2)  w2'2(s2) 

<  C,(l*l  +  1  ♦  I ‘‘  +  IA*I  .  +  l  +  l  ,  ,  ) 

4  LP(S2)  lp(s2>  l2(s3)  w1'2(s3) 

<  C.(|#I  +  >$ln  +  -  +  »*ln,  +  ,*1  1  2  ) 

5  lp(s2)  lp ( s 2>  l2(s3)  l  (s3>  w  '  (S3) 


<  c(l*«  1  2  +  1*1*  2  ). 

W  '  (S„)  w1,2(s,) 


which  proves  ( 4 . 26 ) , 


A  solution  of  equation  (3.1)  obtained  in  section  3  will  be  relevant  to  a 
physical  flow  only  if  we  can  restrict  its  size  by  restricting  some  system 
parameter.  The  previous  lemma  addresses  this  point.  The  following  corollary 
continues  the  development  in  this  direction.  Both  of  these  results  will  be 
used  in  taking  limits  of  periodic  solutions  as  the  period  grows  indefinitely 
large. 


2 

Corollary  4.5.  Let  N=R  if  I{$)  =  R  is  specified  as  in  theorem  3.2 
and  let 

5-m  _ 

I) 

(4.28) 


N  =  6 


(U-X) 


4 ( m- 1 ) 


\a-a 


if  A  is  specified  as  in  theorem  3.4.  Then 

l<M  1  <  C(N  +  Nn) 
C 


and,  for  some  u>  >  0 


(4.29) 
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(4.30) 


>*«  2>u)  <  C'(N  +  nV. 


Here  C  and  C'  depend  on  u  and  on  the  parameters  entering  condition  A;  in 
addition,  C*  and  (o  depend  on  N  through  the  HOlder  exponents  of  f  and 
h. 


Proof.  From  (3.9),  (3.16)  and  lemmas  4.3  and  4.4,  we  have 

l<t>l  <  C(N  +  Nn)  (4.31) 

C(V 

where  N  =  R  or  N  =  6.  according  to  whether  problem  (PC)  or  (PF), 
respectively,  is  being  considered.  The  range  of  $  is  thus  limited  in  terms 
of  R  or  6 .  It  is  assumed  in  condition  A  that  both  f  and  h  are  locally 
HOlder  continuous.  Hence,  for  the  range  of  value  of  (y,t)  which  are 
relevant  for  the  solution  $,  Af-h  is  a  member  of  some  HOlder  class,  say  with 
exponent  id  >  0.  Regularity  theory  for  the  Laplacian  (Gilbarg  and  Trudinger 
1977)  then  yields  (4.30). 


Lemma  4.6.  Let  L  =  -  A  +  q(y)  with  q  continuous  on  0  4  y  <  1. 

2  2 

Suppose  that  the  smallest  eigenvalue  Y^  of  -d  /dy  +q  on  0<y<1 
with  y ( 0 )  =  y ( 1 )  =  0  satisfies  Y  ^  >  0 .  Then  for  u  e  W2,  2((lk)  n  Hk, 

«ul  4  ClLul  ,  (4.32) 

w2,z(flk)  L  (V 

where  C  depends  on  Y 1 • 


Proof .  From  (3.7)  and  (4.25)  it  follows  easily  that 

4  C  _ I  I Lul  .  +  lul  -  _  ) , 


'n  u# 

3  w2,z(s.) 


C  .("Lul  +  «ul 

L2(S.)  ,  w  ,2(S.) 


(4.33) 


where  .  as  in  (4.23),  and  j  =  -k,  -k+1 , . . . ,k- 1 .  It  is  well  known 

3  3  * 1 
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that  for  each  e  >  0  there  is  a  constant  Cg  such  that 

lul  .  ,  <  e»ul  +  C  lul  , 

W1,2(Sj)  W2,2(Sj)  L2^) 


(4.34) 


(cf.  Agmon,  Douglis  and  Nirenberg  1959,  p.  698).  Using  the  partition  of  unity 
{n _j}  and  the  fact  that  u  is  2k  periodic,  (4.33)  and  (4.34)  may  be 
combined  to  obtain 

2  V  2 

lul  -  I  1  n.ul  2  2 

w2'2<nk)  j— k  3  w2'2(fl) 

k-1 

<2 1  m.ui2 

j— k  3  w  '  (^  ) 

k-1  , 

<  l  0.(1 Lul 2  +  elul2  +  C  lul 2  ) 

j—  k  1  L2(sj)  W2,2(Sj)  e  L2(S_j  ) 

<  20,(1  Lul2,  +£lul2,  +  C  lul2  )  . 

1  L2(0k)  w  <V  llV 

The  factors  of  2  appear  since  the  rectangles  S ^  (and  correspondingly  the 
supports  of  the  overlap  in  pairs,  including  S-k  and  Sk-1  if  we 

envision  integration  as  taking  place  on  the  cylinder  obtained  from  the  product 
of  0  4  y  4  1  with  a  circle.  The  choice 


e  *  (40.)  yields 


lul  <  C  (iLul  +  lul  ). 

W2'2(«k)  2  L2(f)k)  L (flk) 


(4.35) 


To  absorb  the  term  lul  we  use  the  hypothesis  on  to  write 

L 


Y1,U“  2 

L  («k, 


<  /{  |Vu|  2  +  qu2} 
=*  / (-Au  +  qu)u 


4  I Lul  .  lul 


obtaining  lul  <  Y^llLul  2 .  This  last  inequality  together  with  (4.35) 

Li  L 

provides  (4.32). 

Next  we  obtain  a  lower  bound  on  —  for  a  solution  (A,$)  of  problem 

L 

P,  assuming  A  <  y.  Recall  that  in  problem  (PF)  we  chose  A  <  y  while  for 
problem  (PC)  it  was  shown  in  lemma  4.1  that  A  <  y  was  a  consequence  of  our 
method  of  obtaining  a  solution. 

Lemma  4.7.  Suppose  f  and  h  satisfy  conditions  AI  and  All  of  section 
1.  Let  (A,$)  e  R  *  with  A  <  y  and  $  >  0,  be  a  solution  of  problem 

(P)(3.1).  Then 

l<J>l  „  >  r  >  0, 

L  (V 

where  r  depends  on  y  -  A ,  but  not  on  k. 

Proof .  Let  5  >  0  be  the  eigenfunction  from  (3.18).  Since 
0  <y  )C  =  /< M-A£  +  hQ5) 

=  / (-A*  +  h0*)£ 

=  /( AfQ4>  +  Af  ^(y,$)  -  h^y^m, 

it  follows  that 

*f..(y,$)  ”  h  (y,$) 

/  <K  [(U~A  ) f 0 (y ) - 1 - j ! - ]  =  0. 

However,  the  term  Af  -  h,  is  o{0)  as  $  ♦  0  and  f  <K  >  0  on  ft  ,  so 

11  o  (  k 

the  integral  cannot  vanish  if  $  is  everywhere  smaller  than  some  r  >  0, 
depending  on  the  gap  y  -  A . 


i 


Corollary  4.8.  For  fixed  R  >  0  in  problem  (PC)  or  for  fixed  X  with 
a""^o  <  X  <  u  in  problem  (PF),  the  solution  $  obtained  in  theorem  3.2  or 
theorem  3.4,  respectively,  has  nontrivial  x  dependence  for  all  sufficiently 
large  k. 


Proof.  In  either  problem  we  have  1$!^  *  S  an<* 


>  r  for 


L  <V 


positive  constants  S  and  r,  independent  of  k.  Suppose  $  were 


independent  of  x.  Then 


0  <  r*  <  1*1  \  <  C  /J  4»2  <  ctk-1l#l  <  Ct) 
L 

which  is  impossible  for  large  k. 


In  fact,  the  symmetrized  solutions  obtained  in  section  3  have  decay  in 
x  over  0  <  x  <  k,  uniformly  in  y,  as  exhibited  in  the  next  lemma. 


Lemma  4.9.  Let  <X,$),  e  be  any  solution  of  problem  (P)  obtained 
from  theorems  3.2  or  3.4.  Then 


♦  (x,y)  <  — *  1/3" - »  if  lxl  < 

I  x  | 1/3 


(4.36) 


where  S  «  l<pl 


Proof .  Let 


6<x)  »  max  $ (x,y) . 

0<y<1 

The  equivalence  of  norms  given  by  (3.9)  combined  with  lemma  4.4  provide  the 


bound 


I  <  C(l$i.  +  l*»“), 

C1  («  ) 
k 


(4.37) 


with  C  independent  of  k.  The  bound  (4.37)  on  the  derivative  of  $  with 
respect  to  y  implies  that  <j>(x,y)  >  6/2  on  a  y  interval  Q  of  length  at 


-41- 


so  for  y  in  Q  and  all  x' 


T 


least  6/C(S+Sn).  But  <P  *  <P  on 
satisfying  -x  <  x'  <  x,  $(x',y)  >  0(x)/2.  It  follows  that 


'♦,22 

L  (8  )  C(S+S  ) 


2  2 

Of  course,  1^1  <  C  S  and  so  a  bound  for  3  results;  viz.  (4.36), 

L  <v 


It  is  possible  to  obtain  faster  decay  in  x.  Bootstrapping  will  furnish 

a  higher  power  of  x  in  the  denominator  of  (4.36).  However,  an  explicit  use 

of  the  Green's  function  for  -A  +  hQ  -  XfQ  will  yield  exponential  decay  for 

♦  ,  as  we  prove  in  the  next  lemma.  For  that  purpose,  let  (Y  »♦  )  n  -  1,2.. 

n  n 

be  the  eigenvalues  and  eigenfunctions  of 


* 


- 1  +  h0(y)i|»  -  XfQ  (y)\|»  =  Y<l», 

dy 


(4.38) 


♦  (0)  =  «|»(1)  -  0, 

with  Y1  <  y  2  <  ^3  ,cnown  (Ince  1927,  p.  270)  that  there  are 

2  2 

constants  CQ  and  C^,  independent  of  n,  such  that  CQn  <  YR  *  C^n  and 

assuming  Y  <  U,  that  Y1  is  at  least  raQ(p  -  Y)  where  mQ  is  the  minimum 

of  fQ(y)  for  0  <  y  <  1.  As  a  normalization  assume  /g'J’2  -  1*  Since 

fl  It-  ♦  I  2  <  (Y  +  Cjjy  for  some  constant  C~,  the  Sobolev  embedding 
Odyn  n  2  '  0  n  2 

inequality  provides  the  crude  estimate,  |  m  <  C,n  for  some  constant  C3 

n  L  J 

independent  of  n.  The  Green's  function  for  -A  +  h^  =  Xf^  in  the  strip  ft 
with  zero  boundary  conditions  on  y  *  0  and  y  =  1  is  easily  computed  using 
separation  of  variables  (cf.  Stakgold  1968,  p.  163).  The  result  is, 


G( x-x' ,  y,  y' )  =  l 


-YV2|x-x'| 

n 


n  2Y 


V2 - VylV»,>- 


(4.39) 


and  one  verifies,  using  the  crude  estimate  on  \|>  ,  that  for  |x-x'|  >  b  >  0 


2- 


(4.40) 


|G(x-x' ,y  y* )  I 


“Y 

<  Ce 


1/2 

1 


|x-x' | 


with  C  depending  on  b. 


Lemma  4.10.  For  k  >  0  let  (A,$),  with  $  e  H^,  be  a  solution  of 
problem  (p)(3.1)  obtained  from  theorem  3.2  or  3.4.  Let  be  the  lowest 

eigenvalue  of  (4.38)  and  suppose  that  for  some  q  >  1 

l^f^y^t)  -  h^ (y,t)  |  <  CQtq 
for  0  <  y  <  1  and  t  >  0.  Then  for  6  <  Y^2/q» 

$(x,y)  <  Ce  ^X^  (4.41) 

and 

|V<fr(x,y)l  <  C’e”^X*  (4.42 

for  |x|  <  k  where  C  and  C*  depend  on  q,  8,  and  l$t  . 

L 

Remark .  For  the  solution  of  problem  (PC)  or  (PF)  we  have  a  bound  on 

l$l  .  Further,  p  -  A  is  specified  or  estimable  and  y  >  m.fp-A)  where 
L  1  0 

mQ  =  inf  fQ.  Hence  $  decays  faster  than  an  exponential  with  a  rate  constant 

1/2 

proportional  to  (p-A)  .  This  dependence  for  the  rate  constant  agrees  with 

that  found  in  the  small-amplitude  theories  of  Ter-Krikorov  (1963)  and  Benjamin 
( 1966). 


Proof.  Since  $  is  positive  and  decreasing  on  0  <  x  <  k  it  is  enough 

to  show  that  the  sequence 

a  =  max  $(i,y),  i  =>  0,1,. ..k, 

0<y<  1 


satisfies 


a.  <  Ce 

l 


-81 


,  i  3  0 , 1 , . . .  k. 


(4.43) 
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< 


We  write  the  equation  (P)(3.1)  as 


(-A  +  hQ  -  XfQ)<j>  =  Xf^y,*)  -  h^y,*). 


and  abbreviate  it  to 


L<j>  =  g  ( $ ) , 


(4.44) 


suppressing  the  dependence  on  X  and  y.  By  hypothesis, 

lg(4»)  I  <  cQ$q. 

Consider  the  cutoff  function  n.  introduced  earlier  in  this  section,  but 

J 

now  extend  it  to  be  2k-periodic  in  x.  Maintaining  the  notation  for  the 


new  function  we  write  equation  (4.44)  as 


i4»  =  n..g(4>)  +  ( 1-n  Jg($). 


(4.45) 


The  operator  L  is  invertible  in  Hk  and  hence  we  can  express  $  as 


^  =  0  +  X  where 


L0  =  n..g(<{>) 


lx  *  (H.)g(^). 


Prom  the  growth  estimate  on  g,  lemma  4.6,  and  the  Sobolev  embedding  theorem 


one  concludes  that 

|0(  j,y)  I  <  c2aj-r 
for  j  =  1,2, ...k.  For  x  one  has  the  estimate 


(4.46) 


IX( j#y) I  <  C  f0  |G( j-x ' ,y  y')l  (1-n  (x' ) )$q(x’ ,y* ) .  (4.47) 

Since  for  each  y,  $  is  even  in  x,  2k-periodic,  and  decreasing  on  the 


interval  0  4  x  <  k,  it  follows  that  for  0  <  j  <  k  and  0  <  x'  <  j, 


<j>(j+x',y)  <  4>(j-x',y). 


(4.48) 


Since  G  is  even  in  its  first  variable,  it  follows  from  (4.48)  that  the 
contribution  to  the  integral  in  (4.47)  taken  over  0  <  x'  <  j  is  at  least  as 
large  as  the  contribution  from  integration  over  j  *  x'  <  2j .  Using  this  last 
observation  and  the  inequality  (4.40)  we  conclude  that 
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ix(  j#y)  I  <  Ct(/j  dx*  +  fpdx'  +  /|x.-j|>jdx,^(1’r,)lGl#qdy’ 
<  0^2  /J  dx'  +  /|x._j|>jdx,)/J(1-n)lGl*qdy* 


jr1 


<  2C  J  /t+1dx'  /q( 1-n)  |G( j-x' ,y#y' ) |$qdy'  + 

*  J  _ A  ^  ^ 


i-0 


(4.49) 


♦  ci  ^|x'-i|>iax‘  (o  lG^-x' *y<y' * 


,  3-1  -r;/2(j-i)  . 

‘  “*  <  io  *  **  *  ~r'/2  L 


Here  we  have  made  essential  use  of  the  vanishing  of  ( 1-n(x' ) )G( j-x' ,y»y’ ) 

-Y  ( 

when  |x'-j|  <  1/4  so  as  to  estimate  it  by  a  suitable  constant  times  e  ' 
for  j-1  <  x*  <  j.  Combining  (4.46)  and  (4.49)  to  estimate  ♦  “  8  +  X  we 
find  that 


1/2 


j-1  -Yl/2(j-i) 


‘j  *  C2aj-1  +  S  \  Jo  * 


a? 


-Y;/2j 


i  Y 1/2 
T1 


l*lq  >  (4.50) 


for  j  -  1,2,... ,k.  Or,  if  aQ  is  redefined  to  be 


1^1 


q  •\1/q 


sup  <$q(0,y)  +  L 


04y<  1 


,V2 


then 


j-i  -Yy2(j-i) 


Sj  <  C  2  aj 


+  c  l  e  1 
-1  3  L  6 


i=0 


(4.51) 


for  j  =  1 , 2, . . . ,k. 
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1  1/2 

Given  8,  0  <  8  <  —  y.  and  B  >  0  we  can  choose  a  positive  integer 
<1  1 

N  depending  on  8  and  B  so  that  for  i  =  0,1,2,...,N, 

8(N-i) 


a^  <  Be 


(4.52) 

To  see  that  there  is  such  an  N  we  first  observe  that  since  $  is  bounded 
and,  by  lemma  4.9,  decays  like 


ai< 


( 1+i) 


1/3' 


for  i  =  0,1,*.. k,  where  E  is  a  positive  constant  depending  on  I$1  The 
function  e^Vtl+i)1^3  is  convex  in  i  and  thus  its  maximum  on  0  <  i  <  N 
occurs  at  i  =  0  or  i  »  N.  Hence  (4.52)  will  hold  provided  N  is  chosen  so 
that 

max( 1 ,  e0N/( 1+N) 1/3 )  <  |  e  . 

E 

We  now  proceed  inductively,  supposing  that  for  some  j  >  N  (4.52)  holds 
for  i  =*  0,1,..., j-1.  From  (4.51), 


*.  <  ♦  c,  Y 

3  2  3  i=0 

1/2 

-  b.® ♦  c  Y  e  ’ 

3  i=0 


) 


8  Y1/2  YV2 

<Be0(N-3)[C2B^V+C3B^1e  1  (1-eq  1  f1]. 


(4.53) 


If  B  is  chosen  so  that  the  contents  of  the  square  bracket  in  (4.53)  add  to 
1,  then  the  desired  inequality  (4.52)  is  valid  for  i  =  j  and  induction 

shows  it  to  be  valid  for  i  =  0,1,. . . ,k.  Thus  (4.43)  and  (4.41)  hold  with 

_  _  8n 

C  =  Be  . 

The  decay  (4.42)  for  follows  readily  from  (4.41)  in  view  of  (4.22), 

(4.27),  and  the  Sobolev  embedding  theorem. 
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5.  THE  EXISTENCE  OF  SOLITARY  WAVES 

In  this  section  we  show  the  existence  of  solitary-wave  solutions  of 
problems  (PC13.3)  and  (PF)(3.4).  These  solutions  will  be  obtained  from  the 
foregoing  periodic  wave  trains  in  the  limit  as  the  period  becomes  unbounded. 


Theorem  5.1.  Suppose  f  and  h  satisfy  condition  A  of  section  3  and 

let  a,  a,  m  be  as  defined  therein;  let  u  be  the  lowest  eigenvalue  of 

(3.17).  Suppose  further  that  m  <  5,  2a  a  <  y  and 

|Xf(y,t)  -  h(y,t) |  <  Ctq  for  some  q  >  1. 

I.  Let  R  >  0  be  given.  For  each  k  >  0,  let  (X^,^)  be  a  solution 

2 

of  problem  (PC)  given  by  theorem  3.2  with  I($  )  *  R  ,  where  I($)  is  given 
by  (3.14).  Then  there  is  an  increasing  sequence  of  half-periods  k(j)  ♦  “, 
j  =  1 » 2, • • • ,  and  a  solution  (X,$)  of  (P)  satisfying 


1) 

4>  > 

0  and 

D 

-e- 

on 

n. 

2) 

I(*) 

% 

Ol 

tc 

It 

3) 

1 4>B 

c 

,2.»  ‘  C<R  *  R” 

)“ 

(cf.  4.30), 

4) 

m 

<  Ce“P 

|x| ;  |V* 

1  < 

Ce"B*X*  (cf  4.41, 

5) 

as 

j  ♦  •» 

Xk(j)  + 

0  <  X  <  y,  and 

6) 

as 

j  ♦  •# 

*k(  j) 

converges  to  $ , 

uniformly  in  C*  on  bounded  subsets  of  (3. 

II.  Let  X  in  (a  ^a,y)  be  given.  For  each  k  >  0  let  (X,$^)  be  a 
solution  of  (PF)  given  by  theorem  3.4.  Then  there  is  a  sequence  k(j)  ♦ 
j  a  1»2,*»»,  as  in  part  I  and  a  solution  (X,$)  of  (P)  satisfying  properties 
1),  4)  and  6)  of  part  I  as  well  as 

2')  H<M  <  ctlim  ■♦■km*km>  < 

W  ,2(0)  j-H»  k(J)  k(;,) 
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and 


3*)  l+l  ,  . .  <  C{6  +  6  )  , 
C2' 


where  6  is  given  by  (4.30). 


Proof .  We  prove  part  I.  The  proof  of  part  II  is  almost  identical. 

Consider  the  sequence  k  =  1,2,...,  restricted  to  ft^.  The  estimate 

(4.30)  for  some  w  >  0  and  the  Arzela-Ascoli  lemma  imply  that  a  subsequence 

2  ( 1 ) 

$H#  $12*  ^13'***'  converges  in  C  (Q^)  to  a  function  $  .  Similarly 

♦  n#  ^12'***'  restricted  to  ft2  contains  a  subsequence  4>  $22'***' 

2  ( 2) 

converging  in  C  (ft to  a  function  $  which,  of  course,  coincides  with 

♦  ^  on  ft^  Continuing  the  process  and  extracting  the  diagonal  sequence 

we  obtain  a  sequence  which  we  denote  3  =  1/2,**»  .  Letting 

♦  be  the  element  of  C  (ft )  defined  to  be  ♦  on  ft  ,  it  is  clear  that  6 ) 

n 

holds.  Without  loss  of  generality  we  can  assume  that  k(j)  is  a  sequence  of 
positive  integers  for  which  k(j)  >  j  and  X^^j  converges  to  a  value  X. 
Since 

W(1  -  C(R) )  <  A  j  <  y ( 1  -  n(R)) 

according  to  lemmas  4.1  and  4.2,  X  has  the  same  bounds,  so  assertion  5) 
holds.  Moreover,  (X  )  is  a  classical  solution  of  the  problem  (P), 


+  h(y,$)  =  Xf(y,$), 
»(x,0)  =  $(x,1)  =  0, 


since  each  pair  (X^j)'  is  a  solution  for  3  =  1/2,...  and  all  terms 

in  the  equation  converge  uniformly  on  bounded  sets. 

Since  $  and  V4  converge  uniformly  on  bounded  sets,  property 

]  ) 

4)  follows  immediately  from  lemma  4.10.  Further  <p  -  <p  on  ft  since,  for 

each  j,  <>,  .  ..  is  even  in  x,  nonnegative,  and  nonincreasing  on 
K  v  J  J 
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0  <  x  <  k(j).  The  estimate  l$l  >  r  >  0  also  follows  from  the  uniform 

L 

convergence  on  in  the  light  of  lemma  4.7.  So  $  >  0  follows  exactly  as 

in  the  proof  of  theorem  3. 2.  Assertion  2)  derives  from  the  convergence  of 
in  C1,  uniformly  on  bounded  sets,  together  with  inequality  4). 
Finally,  the  HOlder  estimates  from  lemma  4.5  are  preserved  in  the  limit, 
yielding  3)  and  completing  the  proof. 


‘  •  ' 


6.  EXAMPLES 


The  application  of  the  foregoing  theory  is  exemplified  in  the  present 
section  for  a  specific  class  of  stratified  fluid  flows.  The  coordinate  system 
and  notation  remain  as  set  forth  in  section  2.  It  is  assumed  that,  in  the 
undisturbed  state  of  the  system,  the  density  p  is  given  by, 

P(y)  =  PqO  “  ay)r,  for  0  <  y  <  1,  (6.1) 

where  p^  and  r  are  positive,  and  0  <  a  <  1.  This  family  of  density 
profiles  is  supplemented  with  a  uniform  primary  velocity, 

U(y)  =  c  ,  for  0  <  y  <  1. 

P 

We  search  for  waves  of  permanent  form  whose  velocity  of  propagation  downstream 
(in  the  direction  of  increasing  x)  is  c.  The  crests  of  the  waves  are 
brought  to  rest  by  considering  the  system  in  a  frame  of  reference  moving  to 
the  right  at  speed  c.  In  the  moving  frame  of  reference,  the  primary  fluid 
velocity  is 

W(y)  =*  U(y)  -  c  =  c^-c=c,  (6.2) 

say.  From  the  specifications  (6.1)  and  (6.2),  the  functions  denoted  by  f 
and  h  in  earlier  sections  may  be  determined  explicitly  in  a  straightforward 
manner  (cf.  Benjamin  1971,  §6.1)  which  is  now  outlined.  The  outcome  of  this 
calculation  is  the  guarantee  that,  subject  to  some  restriction  on  the 
parameters  a  and  r,  the  system  that  has  the  primary  specifications  (6.1) 
and  (6.2)  can  support  both  solitary  waves  and  periodic  waves  of  permanent 
form,  the  latter  being  the  analogue  of  the  classical  cnoidal  waves. 

To  begin  the  computation  of  f  and  h,  note  that  the  pseudo-stream- 
function  associated  to  the  primary  flow  is 


't' (y)  *  /q  p1/2(z)w(z)d2 


C  p 


1/2 


a[1  +  (r/2)  ] 


[1  -  (1  -ay>1+!r/2)]. 

.1/2/ 


It  will  be  convenient  to  let  s  =>  1  +  (r/2)  and  Y  =  c  p^  /as.  With  this 
notation,  we  have 

*(y)  -  Y[1  -  (1  -ay)8].  (6.3) 

Note,  for  later  use,  that 

*yy(y)  «  -s(s-1)a2V(1  -  ay)8”2.  (6.4) 

As  all  variables  need  eventually  to  be  cast  in  terms  of  the  pseudo-stream- 
function,  it  is  convenient  to  express  y  as  a  function  of  't.  The  expression 
is, 

y 

Then  p  may  be  expressed  in  terms  of  f  as 


V(f)  -  4  (1  *  (1  *)1/8]. 

a  T 


(6.5) 


In  particular. 


pm  -  P0d  -  £  »)r/8. 


p*m  «  -  97-  <1  -  7  t> 


(6.6) 


(6.7) 


Similarly,  combining  (6.4)  and  (6.5),  it  follows  that 

Tyy  -  -s(s-1)a2T(1  -  ~  T)(b~2)/b.  (6.8) 

Remember  the  dynamical  condition  satisfied  by  the  flow,  that  the  total  head, 

h  -  p  ♦  gyp  ♦  ~  f*,  (6.9) 

is  conserved  along  each  pseudo-stream- line.  Thus  H  is  a  function  of  f 
alone,  as  pointed  out  already  in.  section  2.  Differentiating  (6.9)  with 
respect  to  V,  and  using  the  hydrostatic  law. 


-gP< 


there  appears 


H*m  =  gyp'd)  +  T  , 

yy 


(6.10) 
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which  is  (2.5)  in  the  special  case  of  a  parallel  flow.  This  completes  the 


calculations  regarding  the  primary  flow. 

Let  4*  denote  the  pseudo-stream-function  for  the  perturbed  flow,  and 

write 

<P(x,y)  =  ¥(y)  +  j<P(x,y) ,  (6.11) 

where  j  is  a  constant,  to  be  fixed  presently.  (It  may  be  either  positve  or 
negative,  corresponding  to  waves  of  depression  or  elevation,  respectively. ) 
According  to  (2.5),  ip  satisfies  the  equation 

-Ai|>  +  H'(»J>)  =  gyp’(40.  (6.12) 

Now,  (6.10),  (6.11)  and  (6.12)  combine  to  yield  the  equation, 

+  f  (*  +  j*)  -  »  (f)  «  -gp*(»  +  j*){Y(»  +  j<t> )  -  Y(t)}, 
satisfied  by  $  (cf.  Benjamin  1971,  eqn.  6.10).  Making  use  of  (6.7)  and 
(6.8),  this  may  be  worked  out  completely  as 

-Aiji  +  h(y,<J>)  =  Xf{y,4>),  (6.13) 

2 

where  A  »  |g  pQr/as  jY It 


f(y,$)  =  s(1  -  ay)S  1sgn(j)f( — ^4 - -) , 


Y  ( 1-ay) 


f(z)  =  (1  -  z)(s'2)/s[1  -  (1-z)1/s], 


(6.14) 


and 


h  ( y  ,4> )  = 


s( s- 1 )Ya' 


( 1  -  ay)  h( 


s-2  rr  — iL 


Y ( 1-ay ) 

_  M  _  M  -  ,(S-2)/s, 


(6.15) 


The  selection  |j|  *  Y  simplifies  the  foregoing  formulae  and  renders  $ 

dimensionless.  If  j  =  1 ,  the  formulae  provide  waves  of  depression,  in  which 

the  streamlines  are  displaced  from  their  undisturbed  configuration  toward  the 

lower  boundary,  whilst  j  =  -Y  corresponds  to  waves  of  elevation,  in  which 

the  streamlines  are  displaced  upward  from  their  undisturbed  postion.  (If 
1/2 

j  =  Y»  then  p  v  ■  -<>  =  -Y$x,  where  v  is  the  vertical  velocity  of  the 

flow.  Considering  the  shape  of  $,  as  established  by  our  theory,  it  is 
concluded  that  the  vertical  velocity  is  negative  to  the  left  of  the  wave's 
crest,  corresponding  to  a  wave  of  depression.  The  opposite  conclusion 
obviously  holds  if  j  **  -Y.)  For  the  sample  primary  flow  described  by  (6.1) 
and  (6.2),  these  two  distinct  cases  correspond  to  different  ranges  of  the 
parameter  r,  and  are  considered  separately  below.  Note  that  if  |j|  *  Y» 
then 

*  * 

c 

and  so  is  proportional  to  the  inverse-square  of  the  velocity  of  the  permanent- 
wave  in  question,  relative  to  the  velocity  of  the  primary  flow. 

It  is  worth  remembering  that  equation  (6.12)  has,  as  an  essential 
ingredient  in  its  derivation,  the  presumption  that  all  the  pseudo-stream-lines 
begin  at  -  *  and  terminate  at  +  •.  It  follows  that  must  take  its 
values  in  the  range  [0,  Y(1)J.  Since  $  vanishes  at  y  *  0  and  y  =  1, 
these  requirements  are  implied  by  the  condition 

*  *  +  j4»y  >  0,  for  (x,y)  e  R  x  (0,1).  (6.16) 

In  the  cases  of  special  interest  here,  this  amounts  to, 

$y(x,y)  >  -sa(1  -  ay)3”1,  (6.17) 

in  case  j  *  Y*  and 

♦  <  sa(1  -  ay)3”1,  (6.18) 

for  j  =  -Y • 
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Waves  of  Depression 

Here  we  concentrate  on  the  case  j  =  Y.  According  to  the  preceding 
discussion,  the  total  head  H  and  the  density  P  are  defined  only  for  values 
of  i|>  lying  in  the  range  [0,  V  ( 1 )  ] .  (In  this  case,  ¥(1)  =  Y[1  -  (1  -  a)S].) 
Hence  the  functions  f  and  h  are  necessarily  as  displayed  in  (6.14)  and 
(6.15)  for  values  of  y  in  [0,1]  and  t  >  0  such  that  ¥(y)  +  Yt  lies  in 
the  interval  [0,  f(1)].  This  implies  exactly  that  the  pair  (y,t)  must 
respect  the  inequality, 

0  <  t  <  (1  -  ay )S  -  (1  -  a)S  =  T(y),  (6.19) 

say.  Note  that  0  <  T(y)  <  t(0)  =V{1)A  <  1,  for  all  y  in  [0,1].  The 
definition  of  f  and  h  for  values  of  (y,t)  outside  the  domain  defined  in 
(6.19)  is  entirely  at  our  disposal.  It  is  technically  convenient,  as  in 
section  3,  to  require  both  f  and  h  to  be  odd  functions  of  the  variable 
t.  So,  only  the  extension  of  the  domain  of  f  and  h  to  encompass  all 
t  >  0  need  be  considered.  For  the  present,  we  content  outselves  with 
extending  f  and  h  to  the  set 

r  =  {(y,t)  :  0  <  y  <  1  and  0  <  t  <  t(0)}  (6.20) 

simply  by  using  the  formulas  in  (6.14)  and  (6.15).  Since  T(0)  <  1,  this 
raises  no  obvious  difficulties.  The  definition  of  f  and  h  for  t  >  T(0) 
will  be  considered  subsequently.  Of  course,  as  discussed  before,  it  is 
required  to  assure,  in  due  course,  that  the  permanent-wave  solutions  obtained 
satisfy  (6.17)  and  so  do  not  depend  on  the  particular  extension  of  f  and 
h  that  is  chosen. 


I 


First,  f  and  h  are  shown  to  satisfy  conditions,  AI,  All  and  AIII,  at 
least  when  restricted  to  the  set  T .  For  this,  we  suppose 


a  <  tt/  [it  +  ( ( s- 1 )  (  2-s ) ) 1/2] 


(6. 21) 


The  restriction  on  a  is  a  technical  one  and  not  a  real  hindrance.  As  r 
varies  between  0  and  2,  s  varies  over  the  interval  [1,2].  Consequently, 
the  constraint  on  a  asks  no  more  than  a  <  .85.  Practical  values  of  a  are 
typically  very  much  smaller.  The  Brunt-V8is&l&  frequencies  N  *  (ag)1/2 
measured  in  the  field  are  of  the  order  .01  to  .02  corresponding  to  an  a 
of  order  10  /cm.  (Xmberger,  Thompson  and  Pandry  1976).  In  the  laboratory, 
Mahony  and  Pritchard  (1981)  measured  N  *  1  or  a  *  .003/cm  for  a  salt 
stratification  near  the  limit  of  salt  solubility. 

It  will  be  convenient  to  represent  f  and  h  as  power  series  in  z  and 

for  that  purpose  to  use 

-v  2  3 

(1-Z)  “  1  +  Vz  +  V(V+1)  jj-  +  v(v+i )  ( v+2)  ~  +  •••  . 

Then,  with 

Pk(5)  -  (5-D5U+1)  •••  (5+k-2),  (6.22) 

the  series  can  be  written 


to 


(1-z) 


1) 


(1-z) 


1) 


00 

-  I 

k=1 


-  vi” fr 


(6.23) 


and 


h(z)  »  -  l 


k=1 


p  (— )  — 

PkV  kl' 


(6.24) 


convergent  for  | z |  <  1 . 

It  is  easily  seen  that,  in  reference  to  hypothesis  AI  of  section  3, 
fQ(y)  *  3tf(y,0)  is  given  by 


fo(y)  *  7^7? 

It  is  strictly  positive  and  infinitely  differentiable  for 
so  is  certainly  locally  Hfilder  continuous  in  that  range. 


(6. 25) 

y  in  [0,1]  and 
In  consequence  of 
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(6.14)  and  (6.23) 


f^y.t)  =  f(y,t)  -  fQ(y)t 


■  ‘V;>  -  V;»  cr  l~~~i 

k=2  (1-ay) 


1.,  1 


. .  .s-l/3-s  1  r 

=  IT?I 


+  J-  [ 


( 1-ay) 


s 


2  31  xS 

s  (1-ay) 


(6.26) 


r  +...} 


=  l  ak(y)fck- 

k=2 


Since  0<t<x(0)<1  it  is  easily  seen  that  f is  locally  HOlder 
continuous  in  T  and  that 

f j(y,t)  <  D(  1  +  t2)  (6.27) 

2  i 

in  T,  for  a  suitably  large  D.  Since  P  (— )  >  0  and  p  (— )  <  0  for 

K  S  K  S 

1  <  s  <  2  and  k=  1,2,3,...,  it  follows  that  ak(y)  >  0  for  all  k.  Thus 
condition  AI  is  seen  to  hold  with 

3-s 

m  =  n  =  2,  a  =  - — ,  and  d  =  D.  (6.28) 

2S 

Regarding  condition  All  relative  to  h,  it  is  easily  seen  from  (6.15) 
and  (6.24)  that 


and  that 


ho(y) 


(s- 1 )  ( s-2)a2 
(1-ay)2 


h.,<y,t)  =  h(y,t)  -  hQ(y)t 


(6.29) 


2  00  P  (— ) 

s(s-1 )a  y  k  s  t  k 

,  ,  .  2-s  k!  ,  ,  s 

(1-ay)  k= 2  (1-ay) 


s(s-1)a 


( 1-ay) 


s+2 


f ( 2-s ) 2  tf  ( 2-s) 2(  2+s ) 

V  2  21  3 

s  s 


3! ( 1-ay ) S 


+ 


(6.30) 

...} 


00 

=  -  l  bk(y)t\ 


where  >  0  for  all  k.  Clearly  hQ  and  h1  are  locally  HBlder 
continuous  in  T.  Since  >  0,  <6.30)  implies  that 

o't2  <  h1 (y,t)  <  ot2  (6.31) 

in  I\  where  0  *  0  and  o'  is  a  suitably  chosen  constant.  By  choosing  a 
possibly  larger  value  of  D  in  ( 6 . 27 ) ,  the  analogous  estimate  for 
I h ^ ( y , t ) |  can  be  assured.  Lastly/ 


3h 

at 


(s-1 ) (s-2)a‘ 
< 1-ay) 2 


[1  - 


-2/s 


(6.32) 


( 1-ay) 


3h  2 

so  the  condition  amounts  to  the  restriction  that  the  quantity  on 
the  right-hand  side  of  (6.3  2)  exceed  -it2  provided  (y#t)  6  T*  In  view  of 
(6.19),  the  latter  is  implied  if,  for  0  <  y  <  1, 


<  ,2. 

(1-ay)  y 

A  short  computation  shows  this  restriction  to  be  equivalent  to  the  supposition 
(6.21)  that 

IT 


a  < 


IT  +  [  ( s- 1 )  (  2-s )  ] 


1/2  * 


Thus  the  condition  All  is  seen  to  hold  on  T. 

Finally,  the  hypothesis  AIII  is  addressed.  Recall  the  definitions  of 
F1  and  ,  that 

f‘1  (y»t)  -2  Jjj  f  1  ( y , p ) dp  and  H^y.t)  =  2  Jj  h  (y,p)dp. 

Since  o=o  in  All,  aa  =0.  Hence  it  is  to  be  verified  that  for  any 

X  >  0,  there  is  a  8  in  (0,1)  such  that 

XF^y.t)  -  H  ^ (y , t )  <  9  [Xf  ^ (y,t)  -  h  1  (y,t) ]  t.  (6.33) 

Rather  than  establish  (6.33)  in  T  we  will  show  that  and  h^  satisfy  a 

slightly  stronger  condition,  which,  in  view  of  the  following  lemma,  will 
guarantee  that  for  suitable  extensions  of  f1  and  h1  outside  I*,  (6.33) 


-57- 


1 


will  hold  for  all  t  >  0.  The  extensions  will  be  such  that  AI  and  All  are 
satisfied  as  well. 

Lemma  6.1.  Let  g  be  a  nonnegative  C1  function  defined  on 
0  <  t  <  tQ  for  some  tQ  >  0.  Suppose  there  is  a  constant  0  with 
2/3  <  0  <  1  such  that 

2g(t)  <  9(g'(t)t  +  g(t))  (6.34) 

1  *v 

on  [0,  tg] .  Then  it  is  possible  to  extend  g  to  a  C  function  g  defined 
for  all  t  >  0  and  satisfying 

G(t)  <  6g(t)t  (6.35) 

where 

G(t)  -2  Jjj  g(p) dp.  (6.36) 

Proof.  If  (6.34)  is  integrated  from  0  to  t,  it  yields 

G(t)  <  0g(t)t,  (6.37) 

for  0  <  t  <  t^,  with  G  defined  in  analogy  with  G.  In  fact  all  one  need 
assume  in  the  proof  is  that  the  inequality  (6.37)  holds  for  0  <  t  <  tQ  and 
that  (6.34)  is  valid  at  the  endpoint  t  =  tg.  Let 

g(t)  =  g(tQ)  +  g' (tp) (t-tg)  +  q(t-tg)2,  (6.38) 

with  q  to  be  chosen.  We  need  check  (6.35)  only  for  t  >  tQ,  given  that 
(6.37)  holds.  Since 

G(t)  -  G(t0)  +  2g(t0)(t-t0)  +  g'(t0)(t-tQ)2  +  j  q(t-tQ)3 

and 

9g(t)t  =  8[g(t0)  +  g'(tQ)(t-t0)  +  q(t-tQ) 2] (t-tQ) 

+  9[g(t0)  +  g'(t0)(t-t0)  +  q(t-t0)2)t0, 

the  desired  inequality  takes  the  form. 
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G(t0)  -  flg(tQ)t0  <  (t-tQ)  t(8-2)g(t0)  +  8g'(t0)t0] 
+  (t-tQ) 2[(8-1)g,(t0)  +  9qt0J 
+  (t-tQ)3[(8  -  |)q] . 

In  view  of  our  assumptions,  the  inequalities  will  hold  if 
q  ■  ( 1-0 )g' (tg)/6t0,  thus  completing  the  proof. 


For  each  y  the  function  g(t)  ■  ^f,|(y»t)  -  h^yft),  defined  for 

0  4  t  4  tg  =  t  ( 0 ) ,  is  a  candidate  to  which  we  may  attempt  to  apply  lemma 

6.1.  If  we  show  g  satisfies  (6.34),  then  it  can  be  extended  to  t  >  tQ 

using  (6.38)  so  that  (6.33)  is  satisfied.  In  order  to  maintain  the  condition 
3h  2 

>  -»  ,  which  is  satisfied  on  0  <  t  <  t(0),  extend  h1  linearly  and  add 
the  quadratic  term  to  f.. 

That  is,  for  t  >  tQ,  define 

£,(t)  -  f,(t0)  ♦  f,(t0)(t-t0)  +  X-\j(t-t0)2 


h^t)  =  h^^)  +  h*(t0)(t-t0). 


(6.39) 


where  the  y  dependence  has  been  suppressed  throughout.  One  can  choose  q 
larger,  if  need  be,  to  guarantee  that  AI  continues  to  hold  for  the  extended 


fr 

1. 

The  inequality  (6.34),  with  9  =  2/3,  holds  for  any  power  t  provided 
k  >  2.  Since  (6.34)  is  linear  in  g,  it  holds  for 


Xf1  -  h  =  l  Uak  +  bk)tk, 

a  series  with  nonnegative  coefficients. 

It  has  been  verified  that  the  functions  f  and  h  arising  from  the 
density  profile  P0(1-ay)r  for  0  <  r  4  2  satisfy  condition  A.  According  to 
theorems  3.2  and  3.4  there  are  x-periodic  solutions  ( A , <#» )  of  equation(3. 1) 
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obtainable  by  specifying  R  in  problem  (PC) (3.3)  or  X  in  problem 
(PF)(3.4).  According  to  corollary  4.5,  condition  (6.17)  will  be  satisfied 
provided  we  restrict  R  to  a  range  0  <  R  <  Rq  or  X  to  a  range 
XQ  <  X  <  p.  with  these  restrictions  in  force, 


with 


c  p 


1/2 


(x,y)  = 


a[1-(r/2)J 
1/2 


[1  - 


(1-ay)1+(r/2)  +  <(>(x,y)]. 


(6.40) 


(*P) 


has  range  in  [0,  f(1)];  is  a  solution  of  Yih's  equation  (2.5);  and  provides, 
for  each  k  >  0,  a  train  of  depression  waves,  2k-periodic  in  x.  Since 
m  <  5  and  0*  2a  ‘a  <  u,  lemma  4.1  implies  that  X  <  y  when  R  is 
specified.  Corollary  4.8  then  guarantees  that  the  x-dependence  is  nontrivial 
for  all  sufficient  large  k. 

For  the  same  ranges  of  R  or  X  there  are  solitary-wave  solutions  of 

Yih's  equation.  These  are  still  described  by  (6.40),  the  pair  (X,$)  now 

2 

being  given  by  theorem  5.1.  Since  IXf^  ~  h^|  *  C^t  ,  the  deviation  of  the 
velocity  fields  of  these  solitary  waves  from  the  velocity  field  of  the  trivial 
flow  decay  exponentially  as  |x|  ♦  “. 

For  r  in  the  range  0  <  r  4  2  (or  1  <  s  4  2),  the  nonlinear 
functions  f  and  h  arising  from  p  *  Pg(1-ay)r  satisfy  condition  A  on  a 
domain  T  which  arose  naturally  in  our  analysis,  and  the  applicability  of  the 
results  of  sections  3-5  was  limited  only  by  having  to  insure  that 
0  <  4  I'M).  For  r  >  2  we  no  longer  deal  exclusively  with  power  series  of 

positive  terms  in  expanding  f  and  -h  and  must  therefore  impose  further 
restrictions  on  the  range  of  $(x,y)  to  apply  the  general  theory.  For  2  <  r 


<4  we  continue  to  use  j  =  Y  so  f  and  h  are  described  by  (6.26)  and 


(6.30)  respectively.  Naturally,  fQ  and  hQ  will  not  change.  To  simplify 
matters,  and  to  exhibit  certain  features  of  the  regime  in  which  2  <  s  <  3, 
attention  is  given  to  the  situation  in  which  0  =  3/4  and  0  <  a  <  . 2,  so 
that  (1-a)s  >  1/2.  Since 

-  |  (k-1) I  <  Pk(5)  <  0,  (6.41) 

for  0  <  5  <  1,  and  since  — —  >  1,  it  follows  that 

1-ay 


. 3-s.  2  "  (2t)k 

f1(y'fc)  >  "  8  i,  "4k“* 

k=3 

From  (6.42)  it  is  easy  to  see  that  the  choice 

a  -2=& 

16s 

2 

more  than  suffices  to  guarantee  f^  >  at  on  the  set 


(6.42) 


(6.43) 


3-s, 


r  IK,< 


(6.44) 

(the  choices  here  are  influenced  by  condition  AIII).  The  remainder  of  AI  is 
verified  as  before.  Again  using  (6.41),  we  see  that  hj,  which  is  in  this 
range  given  by  a  power  series  with  non-negative  coefficients,  satisfies 


0 ‘  hi  ‘  [1r  j3  ^ 

( 1-a)  s  k=3 


] , 


and  hence  for  (y,t)  e  T, 


with  a ’  =  0  and 


;,t2  <  h^(y,t)  <  at2, 


4 ( s-1 ) ( s-2)a' 


(6.45) 


Since  3h/3t  >  0  the  remainder  of  All  is  plainly  satisfied. 


If  condition  (6.34)  is  valid  in  T,  for  g  =  lf^  -  h^,  then  appeal  to 
lemma  6.1  assures  that  f^  and  h^  may  be  extended  to  all  of 
[0,1]  *  R  whilst  still  satisfying  condition  AIII.  As  before,  the  extensions 
described  in  (6.39)  will  serve  to  maintain  conditions  AI  and  All.  Inequality 
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(6.34)  takes  the  form, 


l  Kk+ne-2]  (Xa  (y)  -  b  (y)]t*  >  0, 
k=2  *  K 


(6.46) 


using  the  notation  of  (6.26)  and  (6.30).  The  condition  (6.46)  will  follow 
from  the  more  explicit  inequality, 


(6.47) 


MI3.-3,  t2  -  .  I  tk> 

k=3  4k ( 1-a) 


.  a2{(39-2)  ieiuga  t2  .  ^  l  tk)  .  0. 

K” 3  K\  (“fl  ) 


obtained  using  (6.41)  and  the  restriction  2  <  s  <  3.  With  9  =  3/4, 

(k+1)6-2  <  3k/4  and  the  expression  multiplying  X  in  (6.47)  is  at  least 

2  r1  3-s  .3,2  r  . k- 2, 

1  4  2s  (4*  l  (2t)  1 

k=3 

and  so  more  than  t2(3-s)/16s  for  (y,t)  in  T.  Hence  (6.47)  will  follow  if 

2  * 

.  2  r,  3-s  a  (s-1)  (s-2)  ...  v  3  - k-2.,  „  . 

tlX^  ttT)  (1  +  l  4,2t)  )]  >  °*  (6*48) 

k=3 


In  AIII  we  assume  that 


x  >  o  =  64 ( s-1 )  ( s-2) a 


(6.49) 


and  because  of  this,  a  simple  estimate  shows  (6.47)  to  be  true  in  T  when 


(6.48)  holds.  To  apply  lemma  4.1  we  require  o/a  <  p/2;  i.e. 

64 (s-1)  (s-2) a2  ^  p 
3-s  2 


(6.50) 


This  last  condition  also  guarantees  that  the  interval  (0/a,  p)  is  nonempty 
(cf.  theorem  5.1).  With  the  extra  condition  (6.49)  fulfilled  one  can,  as  in 
the  case  0  <  r  4  2  assert  that  for  2<r<4(2<s<3)  in  (6.1)  there 


I 


are  periodic  and  solitary  waves  of  depression  for  a  suitable  range  of 
"energies"  R  or  supercritical  velocities  c,  derivable  from  (6.40). 

In  verifying  the  condition  A  in  the  range  2<r<4(2<s<3)  we  — ? 

forfeited  any  advantage  of  having  s  near  2  to  illustrate  that  the  range  of 
allowable  amplitudes  vanishes  as  s  approaches  3.  This  is  a  manifestation 
of  a  transition  which  takes  place,  from  waves  of  depression  to  waves  of  * 

elevation.  In  the  limit  of  vanishing  stratification  (a  ♦  0)  the  transition 
occurs  at  r  =  4.  For  positive  a  the  transition  occurs  at  an  r  >  4.  This 
can  be  established  by  relaxing  condition  A  and  using  the  techniques  of  section  ® 

3.  For  exammple,  in  problem  (PF)  one  need  put  conditions  only  on  the 
combination  Xf  -  h.  As  noted  earlier,  condition  A  is  a  compromise  aimed  at 
limiting  the  length  of  the  exposition.  * 


Waves  of  Elevation 

Here  we  are  concerned  with  the  profile  (6.1)  with  an  exponent  r  >  4, 

corresponding  to  s  >  3.  The  terras  fQ  and  hQ  are  unaffected  by  the  choice 

of  j  in  (6.14)  and  (6.15),  but  an  inspection  of  (6.26)  shows  that  the  term 

quadratic  in  t  in  -f^(y,-t)  has  a  positive  coefficient  for  s  >  3, 

dictating  the  choice  j  =  -Y  ♦  To  simplify  writing,  and  so  that  we  can  make 
use  of  the  foregoing  formulas,  we  use  a  superscript  +  to  designate  functions 
associated  with  the  range  s  >  3.  For  example,  f|(y,t)  =  -f  (y,-t). 

The  range  restriction  on  t  now  derives  from 

0U-Tt<),(1), 

and  is  easily  reduced  to 

0  4  t  <  1  -  ( 1-ay)S.  (6.51 ) 

For  the  series  (6.26)  to  be  usable,  t/(1-ay)s  cannot  exceed  1.  Requiring 
-f(-z)  (cf.  6.14)  to  have  a  positive  second  derivative  would  impose  a 
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limitation  of  the  same  order  of  magnitude.  In  combination  with  (6.51)  the 


bound  on  t  requires  that 

1/2  <  ( 1-a) S.  (6.52) 

For  s  =*  3  we  know  that  a  <  .2  will  suffice.  However,  for  large  s  (6.52 
implies  that  a  <  3/s.  Even  with  this  restriction,  we  will  see  in  the  next 
subsection  that  the  exponential  profile  fits  within  the  theory,  treatable 
either  as  a  limit  as  s  grows  and  a  decreases,  or  through  a  direct 
computation. 

In  this  section  the  variable  t  is  most  severely  restricted  by  (6.51) 
when  y  is  near  0  (compare  6.19).  While  the  functions  f^(y,t)  = 

-f^fy^t)  and  h*(y,t)  =  -h1(y,-t)  are  well  defined  in  the  domain  described 
by  (6.20),  it  will  be  necessary  to  restrict  the  range  of  t  in  order  to 
satisfy  condition  A.  The  estimates  here  are  simplified  somewhat  by  an 
alternation  in  the  sign  of  the  terms  in  the  series  expansion  of  f  and  g. 

To  begin,  we  show  that  the  terms  in  the  series  for  f*  alternate  in  sign  and 
decrease  in  size  for  k  >  3.  It  will  suffice,  noting  (6.14)  and  (6.26),  to 
show  that  for  k  >  3  and  0  <  5  *  1/s  <  1/3, 

[Pk(C)  -  Pk(2C)]  (6.53) 

is  non-negative  and  non-increasing  with  increasing  k.  Using  (6.22),  one 
reduces  the  question  to  showing  that 

°  <  [Pk(S)<e+k-1)  -  Pk(25)<  25+k-1)J  <  Pk(C)  -  Pk(2E), 


or  that 


(5-1)5(5+1)  »«»  (S+x-2)  <  2-25 

(26-1)25(25+1)  •••  (  25+k- 2)  2-5 


(6.54) 


A  short  computation  shows  (6.54)  to  hold  for  k  =  3.  For  larger  k  the  left 
side  of  (6.54)  becomes  smaller,  so  that  the  desired  result  holds  for  all 
k  >  3.  It  follows  that 


(6.55) 


ft  > 

1  2s 


< 1-ay) 


s+1 


7-s  t _ 

2  ,  ,  .  ,  .  2s+ 1 ' 

s  3 ! (  1-ay) 


for  0  <  t  <  1  -  (1-a)S.  Since  (1-a)s  >  .5,  f*  >  at2  holds  with 


a  = 


s-3 

4s 


in  the  domain, 


r  =  { (y,  t)  |  0  <  y  <  1,  0  <  t  <  — } 

8(s  -7) 


(6.56) 


Verification  of  AI  in  T  is  completed  as  before.  Again  we  see  that  the  range 
of  t  vanishes  as  s  approaches  3. 


Since  Pk(5)/kl  decreases  with  increasing  k  for  £  =  2/s  <  2/3,  one 
sees  from  (6.30)  that  the  series  for  h*(y,t)  =  -h.|(y,-t)  has  alternating 
signs  and  decreasing  size.  Thus 


a't2  <  h*(y,t)  <  at2 


holds  in  T,  for  a  =  0  and 


a'  =  - 


a  (s-1)  (s-2) 


s( 1-a) 


s+2 


The  t  derivative  of  h+  is  merely  (6.32)  with  t  replaced  by  -t  and  is 
always  positive  for  s  >  3,  so  All  can  be  seen  to  hold. 

To  satisfy  condition  AIII  we  again  arbitrarily  take  0  =  3/4  and  note 
that  t  <  2/3  in  T  from  (6.56).  Then  for  k  >  4,  [(k+1)9-2]tk  is 
decreasing  with  k  and  we  can  exploit  the  alternation  in  the  series  for 
Af^  -  h^  in  satisfying  (6.34)  (compare  (6.45),  (6.46)).  Here  A  >  a/a  =  0 
and  the  analogue  of  (6.46)  will  be  satisfied  if 


,  ,1  s-3  2  s2-7  t3 

4  ?  ^  ^  e 

2s  s  3 !  (  1-a ) 


+  postive  terms) 


a  (s-  1 )  (  3-  2)  3  2  2(  2+s) 

s(  1-a)  4  s 


3!  (  1-a)' 


+  positive  terms)  >  0. 
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Since  (1-a)s  >  .5,  one  easily  checks  that  the  last  inequality  holds  in  T 
The  discussion  surrounding  formula  (6.40)  is  now  relevant,  though  with  $ 
replaced  by  -$  in  that  formula,  corresponding  to  waves  of  elevation. 


The  Exponential  Profile 


xf  the  density  in  the  undisturbed  state  is  taken  to  be 

P(y)  =  PQe  ^y 

then  the  pseudo-stream-function  corresponding  to  (6.3)  is 


ap01/2c  -  1  By 

’•'(y)  =  — « —  (l  -  e  ). 


(6.57) 


(6.58) 


A  computation  similar  to  that  carried  out  for  the  density  (6.1)  shows  that  in 
this  case  the  basic  equation  (6.13)  (or  (2.1))  has 

—  —  By  -  By  -  By 

f+(y,4>)  -  e  2  (l  +  e2  <f>)  log(  1  +  e2  $),  (6.59) 


h+(y,«fr)  «  j-  4>, 


\  =  g(J/c  , 


(6.60) 


(6.61) 


with  j  =  -y  =  “2Pq  c/S.  One  can  now  proceed  to  verify  condition  A  as 
before.  However,  we  have  already  done  this,  for  (6.57)  is  obtainable  as  a 
limit 

P.e'0y  =  lim  p  (1  -  -  y)r 

0  r-H»  0  r 

with  a  *  B/r  in  (6.1).  In  terms  of  s,  with  a  =  0/(2s-2),  (1-ay)s 


approaches  e 


-  -  0y 

2  y  ,  + 


and  f  approaches  1  as  s  +  +  ®.  The  expression 


s(Pk(2/s)  -  P  (1/s)]  in  (6.26)  approaches  the  derivative  of  Pk  at  0, 


i 


4 


i.e.  (k-2)l,  and  the  whole  series  in  (6.26)  approaches  a  series  for 
f(y,t)  from  which  one  obtains  a  series  for  f*(y,t)  =  -f^yj-t).  Of  course, 
this  series  can  also  be  calculated  directly  using  (6.59). 


From  (6.29)  it  is  apparent  that 


+  8 
h0(y)  = 


while  from  (6.30),  h^  =  0,  since  each  coefficient  is  Q(1/s),  as  s  +  +  °». 
Previously  we  had  the  restriction  a  <  3/s  for  large  s  so  that  now  it  is 
required  that  0  =  2as  <  .6.  The  restriction  on  t  from  (6.56)  allows  us  to 
take  0  <  t  <  3/8  and  from  (6.55),  a  =  1/4  will  suffice  on  that  interval. 
The  remaining  points  of  condition  A  carry  over  in  the  limit  and  we  therefore 
obtain  periodic  and  solitary  waves  of  elevation  with 


where 


2Pq/2  c  -  1  3y 

4>(x,y)  =  — ^ -  [1  -  e 


c  =  (g0/\)V2. 


<Mx,y)]  , 


(6.63) 
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7.  CONCLUSION 


•  ^ 


The  hypotheses  AI,  All  and  AIII  introduced  in  section  3  and  the 
additional  restrictions  in  lemmas  4.1  and  4.10  effectively  delineate  a  class 
of  quiescent  states  of  a  heterogeneous  fluid  confined  between  two  rigid 
horizontal  plane  boundaries.  Wave  motion  superimposed  upon  these  quiescent 
states  has  been  contemplated.  Provided  that  the  effects  of  dissipation, 
diffusion  of  the  heterogeneity,  and  compressibility  may  be  ignored,  on  time 
scales  sufficient  for  the  passage  of  the  waves  in  question,  this  class  of 
heterogeneous  systems  has  been  demonstrated  to  support  two-dimensional  steady- 
wave  motion  corresponding  generally  to  surface  cnoidal  and  solitary  waves.  It 
deserves  emphasis  that  these  waves,  which  are  solutions  of  the  full  Euler 
equations,  need  not  have  small  amplitude. 

In  this  section  we  discuss  some  implications  and  drawbacks  of  the  present 
theory,  and  point  to  directions  for  future  investigation. 

Perhaps  the  most  obvious  drawback  of  our  analysis  is  the  absence  of  a 
practical  set  of  sufficient  conditions  that  a  given  system  fall  within  the 
confines  of  the  class  covered  by  the  theory.  Such  conditions  would  ideally  be 
e^qpressed  directly  in  terms  of  the  undisturbed  density  stratification  p  and 
the  undisturbed  fluid  velocity  U.  As  the  examples  studied  in  section  6  amply 
demonstrate,  our  theory  does  have  utility  in  a  considerable  range  of 
situations.  Nevertheless,  checking  the  hypotheses  of  the  theory  can  be 
somewhat  tedious. 

It  should  also  be  admitted  that  variation  of  density  that  features  a 

relatively  rapid  change  in  the  middle  of  the  vertical  extent  of  the  flow 

domain,  say  joining  two  homogeneous  layers  of  fluid  of  constant  density 

P„  and  p„  (with  p,„  <  p  )  situated  at  the  top  and  bottom  boundary, 

T  B  T  B 

respectively,  have  not  been  treated  by  our  theory.  Such  density  variations 


•  * 


1 


•  ' 
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are  addressed  in  the  complementary  paper  of  Turner  (1981).  The  latter 
reference  even  includes  results  pertaining  to  the  much-studied  two-layer 
problem  (a  layer  of  light  fluid  resting  upon  a  layer  of  heavier  fluid)  which 
has  a  discontinuous  undisturbed  density  profile.  Similarly,  Ter-Krikorov 
(1963)  was  able  to  treat  systems  having  an  arbitrary  smooth  undisturbed 
density  profile.  Ter-Krikorov' s  techniques  limit  his  results  to  small 
amplitude  solitary  and  cnoidal  waves.  While  the  variational  techniques  used 
by  Turner  allow  for  finite-amplitude  waves,  the  required  estimates  for  the 
quasilinear  elliptic  equation  arising  in  his  analysis  effectively  limit  the 
applicability  to  small  amplitudes.  We  faced  a  similar  situation  here  in 
having  to  limit  the  ranges  of  energies  or  speeds.  However,  the  linear 
elliptic  estimates  intervening  in  section  4  are  not  such  as  to  vitiate  our 
claim  to  have  found  "finite-amplitude  waves".  In  none  of  the  studies  under 
discussion  is  the  entire  collection  of  allowable  solutions  exhausted.  In  the 
case  of  surface  waves  on  fluid  of  constant  density  one  has  a  more  satisfactory 
global  picture.  Amick  and  Toland  (1980)  have  shown  that  there  is  a  connected 
set  of  solutions  joining  a  "wave  of  greatest  height"  to  a  trivial  flow. 

Another  possible  objection  to  the  present  theory  is  the  fixed  upper 
boundary  that  intrudes  in  our  formulation.  This  simplification  has  a  long 
history  (cf.  Keulegan  1953,  Long  1956,  Benjamin  1966,  1971).  Despite  these 
precedents,  a  fixed  upper  boundary  certainly  seems  somewhat  artificial  when 
one  thinks  of  waves  in  the  ocean  or  atmosphere.  Peters  and  Stoker  (1960) 
obtained  a  formalism  for  internal  solitary  waves  which  allowed  a  free  upper 
surface.  Benjamin  (1966,  §4)  also  gives  theoretical  results,  some  of  which 
relate  to  the  situation  where  the  upper  surface  is  left  free.  He  notes  that 
the  resultant  change  in  boundary  conditions  that  is  implied  in  passing  from  a 
fixed  to  a  free  upper  surface  can  fundamentally  alter  the  character  of  the 
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flow.  Especially  considering  this  latter  point,  a  rigorous  theory  that  takes 
account  of  a  free  upper  surface,  or  which  is  set  in  a  semi-infinite  domain 
where  the  upper  boundary  has  receded  to  infinity,  would  seen  a  worthwhile 
undertaking. 

Another  aspect  of  internal-wave  theory  that  deserves  attention  is  the 
relationship  of  solutions  of  the  Euler  equations  to  solutions  of  the  various 
model  equations  that  formally  apply  to  small-amplitude,  long-wavelength 
disturbances.  Because  of  the  extra  length  scale,  the  vertical  width  of  the 
stratified  layer,  a  number  of  essentially  different  models  arise,  a  situation 
at  variance  with  the  usual  theory  of  small-amplitude  long  suface  waves. 
Depending  on  the  relationship  between  the  varius  parameters  that  intervene  in 
the  regime  at  hand,  we  may  obtain  the  Korteweg-de  Vries  equation  (Benney, 
1966),  Benjamin's  (1967)  equation,  or  the  equation  derived  by  Joseph  (1977). 
(The  latter  equation  reduces  to  the  former  two  in  certain  limiting  cases.)  It 
is  not  intended  to  initiate  a  detailed  discussion  of  these  equations. 

However,  it  is  interesting  to  recall  that  the  solitary-wave  solutions  of  these 
models  play  a  distinguished  role  in  the  evolution  of  general  classes  of 
initial  data,  in  that  the  initial  data  resolves  itself  into  a  sequence  of 
solitary  waves  followed  by  a  decaying  dispersive  tail  (cf.  Miura  1976).  There 
are  not  results  available  now  providing  rigorous  comparisons  of  solutions  of 
these  model  equations  with  solutions  of  the  Euler  equations.  Notwithstanding, 
such  model  equations  have  been  observed  to  give  quite  acceptable  predictions 
(cf.  Bona,  Pritchard  and  Scott  1981,  and  the  references  contained  therein),  so 
it  is  not  unreasonable  to  conjecture  that  solitary  internal  waves  do  play  a 
somewhat  special  role  in  the  long-term  evolution  of  certain  classes  of 
disturbances.  (Some  additional  evidence  in  favour  of  this  conjecture  will  be 
introduced  presently.)  If  valid,  this  proposition  certainly  heightens 
interest  in  solitary-wave  solutions  of  the  Euler  equations. 


The  waves  whose  existence  is  demonstrated  in  section  5  deserve,  at  least 


provisionally,  the  appellation  "solitary  wave".  They  accord  with  Scott 
Russell's  (1845)  original  conception,  a  single-crested  wave  which  is  sensibly 
localized  in  space  and  which  propagates  at  a  constant  velocity  and  without 
change  of  form.  The  surface  water  waves  regarded  in  the  field  and  in  the 
laboratory  by  Scott  Russell  are  remarkably  stable.  To  fully  justify  the 
designation  "solitary  wave"  for  the  waves  dealt  with  herein,  a  satisfactory 
theory  of  stability  pertaining  to  these  waves  must  be  forthcoming.  Alas,  at 
this  juncture  even  a  uniqueness  result  may  only  be  surmised.  Therefore,  this 
issue  is  left  aside  in  the  present  work,  though  its  resolution  would  certainly 
be  of  considerable  interest.  The  solitary-wave  solutions  of  the  Korteweg-de 
Vries  equation  and  of  Benjamin's  equation  are  known  to  be  stable  (cf.  Benjamin 
1972,  Bona  1975  and  Bennett  jst.  _al.  1981).  Insofar  as  these  model  evolution 
equations  reflect  the  Euler  equations,  the  last-mentioned  results  may  be 
construed  as  evidence  in  favour  of  the  stability  of  the  waveforms  dealt  with 
in  section  5. 

Experimental  and  field  observations  point  even  more  convincingly  to  the 
stability  of  the  internal  waves  addressed  earlier.  We  point  to  the  laboratory 
studies  of  Davis  and  Acrivos  (1967),  Walker  (1972)  and  Maxworthy  (1979, 

1980).  On  a  grander  scale,  the  recent  observations  of  internal  waves  off  the 
coasts  of  the  continental  U.S.A.  (cf.  Apel  ^t.  al.  1975),  in  the  Andaman  sea 
(Apel  1979,  Osborn  and  Burch  1980),  and  elsewhere,  present  themselves. 

Internal  waves  with  peak  to  trough  amplitudes  as  large  as  100m  are  represented 
in  these  rather  spectacular  reports  from  the  field.  Some  of  these  waves  have 
been  interpreted  as  internal  solitary  waves,  and  indeed  the  evidence  for  this 
presupposition  is  substantial,  if  not  compelling.  From  the  view  of  the 
present  discussion,  these  observations  tend  especially  to  reinforce  two 


'  • 
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points.  The  first  is  the  exceptionally  stable  nature  of  internal  solitary 
waves.  The  second  is  the  distinguished  role  played  by  the  solitary  wave  in 
the  long-term  evolution  of  certain  types  of  disturbances,  as  would  be 
predicted  by  the  various  model  equations.  As  the  waves  in  question  probably 
obey  the  Euler  equations  more  closely  than  the  simplified  model  equations, 

there  is  an  inference  that  the  Euler  equations  may  exhibit  the  forementioned 

* 

phenomenon  of  resolving  perturbations  into  a  sequence  of  solitary  waves  . 

Finally,  it  deserves  emphasising  that  the  periodic,  or  cnoidal  waves 
whose  existence  has  been  assured  in  theorems  3.2  and  3.4  do  not  satisfy  the 
boundary  conditions  under  which  Yih's  equation  (2,5)  was  derived.  In 
consequence,  while  they  may  be  used  to  define  solutions  of  the  Euler 
equations,  no  direct  physical  significance  can  be  imputed  to  these 
solutions.  Nonetheless,  it  has  been  argued  by  Benjamin  (1966,  §3)  that  such 
periodic  wavetrains  may  be  generated  in  the  lee  of  an  obstacle  introduced  into 
a  primary  flow,  or  on  the  back  of  an  internal  undular  bore.  In  both  these 
cases,  a  change  in  energy  is  needed  in  the  generation  region,  near  the 
obstacle  or  near  the  front  of  the  bore,  respectively,  in  order  that  such  waves 
may  arise.  In  consequence,  this  sort  of  situation  is  beyond  the  scope  of  the 
present  study,  both  in  regard  to  predicting  which  wavetrain  would  be 
manifested  in  a  given  configuration,  and  in  regard  to  the  possibility  of  an 
exact  analyis.  In  any  case,  the  periodic  wavetrains  present  some  interest,  if 
for  no  other  reason  than  their  approach  to  the  solitary  wave,  as  their  period 
grows  indefinitely  large. 

*  '  ’  - 

We  are  not  suggesting  that  the  Euler  equations  are  completely  integrable. 

The  recent  work  of  Benjamin  and  Olver  (1981)  makes  this  an  unlikely 
possibility.  However,  it  appears  that  wave  equations  may  have  the  property  of 
resolving  initial  data  into  solitary  waves  without  being  completely  integrable 
(cf.  Bona  1981). 
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Appendix 


•  ' 


In  this  section  we  set  forth  some  results,  used  in  section  3,  regarding 
symmetrization  (or  rearrangement)  of  functions* 

Let  T  be  a  triangulation  of  the  rectangle  M  {(x,y)  |  -k  <  x  <  k, 

0  4  y  <  l},  i.e.,  a  collection  of  closed  triangles  contained  in 

having  union  equal  to  and  having  the  property  that  any  pair  intersect  in 

a  common  edge  or  vertex,  if  at  all.  Let  u  be  a  function  defined  on  We 

say  u  is  piecewise  linear  (PL)  if  it  is  continuous  in  and  if  there  is  a 

triangulation  T  of  as  described  above,  such  that  u  is  linear  in  x 

and  y  on  each  triangle  of  T.  Starting  from  the  fact  that  C1^)  is  dense 
in  W  '  (tl^)#  and  approximating  a  C  function  by  a  PL  function  in  the 
W' z  norm,  one  shows  that  PL  functions  are  dense  in  W  '  (Q^).  We  denote  by 

PL+  the  collection  of  u  6  PL  for  which  ux  ?  0  a.e.  Suppose  s  e  PL  has 

sx  +  0  a.e.  By  considering  u(x,y)  +  es(x,y),  for  suitable  small  e,  one 
sees  that  PL+  is  dense  in  W^'^(O^).  All  we  require  for  our  variational 
principles  is  a  dense  class.  At  the  same  time,  the  symmetrization  which  we 

next  define  is  very  easily  analyzed  in  the  class  PL*.  One  can  extend  the 

1  2 

results  on  symmetrization  to  all  of  W  '  as  was  done  by  Fraenkel  and  Berger 
(1974).  However,  the  weak  limits  involved  can  equally  well  be  postponed  until 
working  with  a  variational  principle  as  we  do  in  this  paper. 

Suppose  u  is  in  PL+  on  For  each  fixed  y  with  0  <  y  <  1,  let 

y(a,y)  =  m{x  |  u(x,y)  >  a}, 

where  m  is  1-dimensional  Lebesgue  measure  and  a  lies  between  the  minimum 
and  maximum  of  u.  The  measure  in  question  is  merely  the  sum  of  the  lengths 
of  a  finite  set  of  intervals  and  since  ux  ^  0  a.e. ,  u  will  be  continous 
in  a  and  y.  For  each  y  the  value  of  y  strictly  decreases  as  a 
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increases .  For  a 


max  U/  (i  ■  0,  while  if  a  =  min  u,  u  =  2k.  We  seek  a 


function  u  which  is  even  in  x,  decreases  in  x  for  0  <  x  <  k,  and 
satisfies 

a* 

m{x  |  u(x,y)  >  a}  =  u(a,y). 

If,  for  a  given  y,  a(u,y)  denotes  the  function  inverse  to  W(a,y),  thi 
the  even  function  u(x,y)  defined  by 

u(x,y)  =  a(x,y), 

for  0  <  x  <  k,  clearly  meets  our  criteria.  Note  that  if  v  =  v(x,y)  is 

A/ 

another  PL  function,  then  u  <  v  implies  u  <  v. 

Lemma  1.  Let  u  and  w  be  in  PL+.  Then 

<v 

lu-wl  m  <  iu-wl  m  . 

L  (flk)  L  (ftk) 


Proof.  Let  lu-wl  _  *  6.  Then  w  <  u  +  6  so  w  <  u+6.  If  u(x)  =  < 

”  —■  OB 

L  _  _ 

then  2x  =  m{u  >  a}  =»  m{u  +  6  >  a  +  6},  so  (u+6)(x)  *  a+6  =  u(x)  +  6.  T1 

“■v  »v  <v  <v 

w  <  u+6.  similarly,  u  <  w  +  6,  so  lu-wl  —  <  5. 

L 


For  u  e  PL+  the  absolute  value  |u|  is  in  PL+.  Define  (compare 
definition  2.1) 


v  -  |  u  | 


and 


u  *  v. 

A 

We  call  u  the  symmetrlzatlon  of  u,  after  Steiner  (cf.  Polya  and  Szego 
1951,  note  A).  It  follows  immediately  from  lemma  1  that 

A  A 

lu-zl  _  <  lu-zl  (Ap.  1) 

00  00  * 

L  (£)k)  L  (ftR) 


for  PL+  functions  u  and  z  defined  on  ft.  . 


Our  use  of  symmetrization  hinges  on  the  behavior  of  certain  integrals  of 
symmetrized  functions.  First,  suppose  G(y,u)  is  continuous  in  y  and  u, 
for  0  <  y  <  1  and  any  real  u.  For  u  e  PL+  an  elementary  argument  shows 
that 

.  u 

L  G(y,u(x,y) )dxdy  -  -J  /  “*  G(y,a)dy (a,y)dy, 
k  min 

and  hence  that 

Jfi  G(y,u) dxdy  -  /_  G(y,u)dxdy 
k  k 

A 

(cf.  Polya  and  Szego  1951,  p<  185).  Since  the  function  u  has  the  same 
measure  distribution  y  as  |ul,  it  follows  that 

A 

Ja  G(y, |u| )dxdy  =  G(y,u)dxdy. 

Mk  k 

If  G(y,|u|)  *  G(y,u) ,  i.e.  if  G  is  even  in  u,  then 

A 

/o  G(y,u)dxdy  =  Jn  G(y,u)dxdy.  (Ap.  2) 

*  k  k 

Another  relation  we  need  for  a  function  u  in  PL+,  which  is  2k 

periodice  in  x,  is  the  inequality 

/n  |Vu| 2dxdy  >  |Vu|2dxdy,  (Ap.  3) 

k  k 

2  2 

which  follows  from  (Polya  and  Szego  1951,  p.  186)  since  |Vu|  =  |V|u| | 
a.e.  It  should  be  noted  that  (Ap.  3)  will  not  hold  for  an  arbitrary  PL+ 
function  on  For  the  cited  proof  to  apply  it  must  be  that  each  line  y  = 

constant,  u  =  constant  which  intersects  the  graph  of  u,  intersects  it  in  at 
least  two  points.  A  function  which  is  2k  periodic  in  x  has  this  property 
on  and  since  we  consider  only  such  functions,  (Ap.  3)  is  certainly  at 


our  disposal  throughout. 
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